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In two recent papers [J, [3 we introduced "dynamical dark matter" (DDM), a new framework 
for dark-matter physics in which the requirement of stability is replaced by a delicate balancing 
between lifetimes and cosmological abundances across a vast ensemble of individual dark-matter 
components whose collective behavior transcends that normally associated with traditional dark- 
matter candidates. We also presented an explicit model involving axions in large extra spacetime 
dimensions, and demonstrated that this model has all of the features necessary to constitute a viable 
realization of the general DDM framework. In this paper, we complete our study by performing 
a general analysis of all phenomenological constraints which are relevant to this bulk-axion DDM 
model. Although the analysis in this paper is primarily aimed at our specific DDM model, many of 
our findings have important implications for bulk axion theories in general. Our analysis can also 
serve as a prototype for phenomenological studies of theories in which there exist large numbers of 
interacting and decaying particles. 
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I. INTRODUCTION 

Dynamical dark matter (DDM) [l|, |2] is a new framework for dark-matter physics in which the requirement of 
stability is replaced by a delicate balancing between lifetimes and cosmological abundances across a vast ensemble of 
individual dark-matter components. Due to the range of lifetimes and abundances of these components, their collective 
behavior transcends that normally associated with traditional dark-matter candidates. In particular, quantities such 
as the total dark-matter relic abundance, the proportional composition of the ensemble in terms of its constituents, and 
the effective equation of state for the ensemble possess a non-trivial time dependence beyond that associated with the 
expansion of the universe. Indeed, from this perspective, DDM may be viewed as the most general possible framework 
for dark-matter physics, and traditional dark-matter models are merely a limiting case of the DDM framework in 
which the states which compose the dark sector are taken to be relatively few in number and therefore stable. 

In Ref. [l|, we laid out the general theoretical features of the DDM framework. By contrast, in Ref. 2], we presented 
an explicit realization of the DDM framework: a model in which the particles which constitute the dark-matter 
ensemble are the KK excitations of an axion-likc field propagating in the bulk of large extra spacetime dimensions. 
We demonstrated that this model has all of the features necessary to constitute a viable realization of the general DDM 
framework. In this paper, we complete our study by performing a general analysis of all phenomenological constraints 
which are relevant to this bulk-axion DDM model. Although the analysis in this paper is primarily aimed at our 
specific DDM model, many of our findings have important implications for theories involving large extra dimensions 
in general. Furthermore, our analysis can also serve as a prototype for phenomenological studies of theories in which 
there exist large numbers of interacting and decaying particles. 

It is important to emphasize why a general analysis of this sort is necessary, given the existence of numerous 
prior studies of the phenomenological and cosmological constraints on axions and axion-like fields, unstable relics, 
and the physical properties of miscellaneous dark-matter candidates. As discussed in Refs. [l|, |2[, such studies are 
typically applicable to dark sectors involving one or only a few fields, and are typically quoted in terms of limits on 
the mass, decay width, or couplings of any individual such field. It is not at all obvious how such bounds apply to 
a DDM ensemble — a dark-matter candidate which is not characterized by a well-defined single mass, decay width, 
or set of couplings. For example, constraints on a cosmological population of unstable particles derived from big- 
bang-nucleosynthesis (BBN) considerations or bounds on distortions in the cosmic microwave background (CMB) 
are generally derived under the assumption that such a population comprises but a single particle species with a 
well-defined lifetime and branching fractions. Such constraints are not directly applicable to a DDM ensemble (in 
which lifetimes are balanced against abundances), and must therefore be reexamined in this new context. 

In this paper, we shall develop methods for dealing with these issues and for properly characterizing the constraints 
on models in which the dark-matter candidate is an ensemble of states rather than a single particle. As we shall find, 
the presence of non-trivial mixings among the KK excitations in our DDM model gives rise to a number of surprising 
effects which arc ultimately critical for its phenomenological viability. One of these is a so-called "decoherence" 
phenomenon 1-3] which helps to explain how the dark matter in this model remains largely invisible to detection. 
Another is a suppression, induced by this mixing, of the couplings between the lighter particles in the dark-matter 
ensemble and the fields of the Standard Model (SM). As we shall see, these effects assert themselves in a variety of 
phenomenological contexts and play a crucial role in loosening a battery of constraints which would otherwise prove 
extremely severe. 

This paper is organized as follows. In Sect. |lTl we briefly summarize the physics of axions in extra dimensions 
and review the notational conventions established in Ref. [4I, which we once again adopt in this work. In Sect. IIIIl 
we examine a number of processes, both thermal and non-thermal in nature, which contribute to the generation of 
a cosmological population of axions. We calculate the rates associated with these processes and assess the relative 
importance of the associated production mechanisms within different regions of model-parameter space. In Sect. II V[ we 
then discuss the phenomenological, astrophysical, and cosmological constraints relevant for bulk-axion DDM models 
and assess how the parameter space of our model is bounded by each of these constraints. In Sect. |Vl we summarize 
the collective consequences of these constraints on the parameter space of our bulk-axion DDM model. Finally, in 
Sect. lVIi we discuss the implications of our results for future research. 



II. GENERALIZED AXIONS IN EXTRA DIMENSIONS: A REVIEW 

In this section, we provide a brief review of the physics of generahzed axions in extra dimensions. (More detailed 
reviews can be found in Refs. Q, 3-) ^Y "generahzed axion," we mean any pseudo-Nambu-Goldstone boson which 
receives its mass from instanton effects related to a non-Abelian gauge group G which confines at some scale Ag. 
Note that the ordinary QCD axion [J, Q is a special case of this, in which G is identified with S'[/(3) color and Kg 
is identified with Aqcd ~ 250 MeV. However, in this paper, we shall leave these scales arbitrary in order to give our 
analysis a wider range of applicability. We will also assume the existence of a global Abelian symmetry U{l)x which 
plays the role played by the Peccei-Quinn symmetry C/(1)pq in the specific case of a QCD axion. 

Our goal in this paper is to study the phenomenological constraints that arise when a generalized axion is allowed 
to propagate in the bulk [3| of a theory with extra spacetime dimensions iS, lZi]. In particular, we consider the case 
in which the axion propagates in a single, large, flat extra dimension compactified on a Sxj'S.i orbifold of radius R. 
The fields of the SM are assumed to be restricted to a brane located at xs = 0. We also assume that the additional 
non-Abelian gauge group G is restricted to the brane at x^ = 0. At temperatures T ^ Aq, the effective action for a 
bulk axion in five dimensions can be written in the form 
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-duad a + Six^) (£branG + -Cint) 



(2.1) 



where 'a' denotes our five-dimensional axion field, ^Cbrano contains the terms involving the brane fields alone, and Ci^t 
contains the interaction terms coupling the brane-localized fields to the five-dimensional axion. The second of these 
terms is given by 
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where F^^, G" , and t/'Jj^ are the field strengths respectively associated with the C/(I)em, SU{Z) color, and G gauge 
groups; -F)ji/, G" , and t/° are their respective duals; e, gs^ and go are the respective coupling constants for these 
groups; fx is the fundamental five-dimensional scale associated with the breaking of the U{l)x symmetry; c^, Cg, 
and Ci are coefficients which respectively parametrize the coupling strength of the five-dimensional axion field to the 
photon, gluon, and fermion fields of the SM; and ^ is an 0(1) coefficient which depends on the specifics of the axion 
model in question. Note that in Eq. (12. 2[) . we have displayed terms involving only the light fields of the SM (i.e., the 
photon, gluon, and light fermion fields), as couplings to the heavier SM fields will not play a significant role in our 
phenomenological analysis. 

The five-dimensional axion field can be represented as a tower of four-dimensional KK excitations via the decom- 
position 
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where the factor 
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(2.4) 



ensures that the kinetic term for each mode is canonically normalized. Substituting this expression into Eq. (J2.2I) and 
integrating over x^,, we obtain 
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where the axion potential is given by 



V{a) = Y. 
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1 n^ 2 
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(2.6) 



and where the quantity fx, defined by the relation 



/i ^ 27ri?/i , 



(2.7) 



represents the effective four-dimensional C/(l)js:-breaking scale. Note that each mode in the KK tower couples to the 
SM fields with a strength inversely proportional to fx- Also note that at these scales, the axion mass-squared matrix 
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(a) 



is purely diagonal. 

At scales T < Aq, an additional contribution to the effective axion potential arises due to instanton effects. In this 
regime, the potential is modified to 






2i?2«"+32^2^G 



1 - cos I ^^ ^ r^an + Qg 
\jx n=0 J 



(2.9) 



where Qq is the analogue of the QCD theta-parameter 0. This results in a modification of the axion mass-squared 
matrix to 
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where Mc = l/R is the compactification scale. This matrix above takes the form [3| 
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The eigenvalues A^ of this mass-squared matrix are the solutions to the transcendental equation 

cot = A^ . 

y \mxyj 

The corresponding normalized mass eigenstates ax are related to the KK-number eigenstates a„ via 



a\ 



y^ Uxnan = ^ 



where A = X/nix- The dimensionless quantity A\ is given by 
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and obeys the sum rules [3| 
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For T <Si fx, rewriting Eq. (|2.5[) in terms of the a\ and expanding the axion potential given in Eq. (|2.9I) out to 
0(a^//^) yields the effective action 

J L A ^ S2n^fx 32n^Jx 
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Of course, the interaction term between the a\ and the gluon field is only a useful description of the physics at 
temperatures above the quark-hadron phase transition at T ~ Aqcd- At temperatures below this threshold, this 
interaction term gives rise to an effective Lagrangian containing interactions between the a\ and various hadrons, 
including the proton p, the neutron n, and the charged and neutral pions tt and tt^ . This Lagrangian takes the form 
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where the coefficients Ca-n, Cam eic, depend on the details of the theory. For example, for a "hadronic" QCD axion Q 
(i.e., a QCD axion which does not couple directly to the SM quarks), the coefficients Cap and Cam which determine 
the strength of the axion-nucleon-nuclcon interactions, are 

Cap = 0.24 (yT^) + 0.15 (^) + 0.02 , Can = 0.24 [r^}) + ^-^^ (tTt) + ^'^^ ' ^^'^^^ 

where z = mu/rnd ~ 0.56 is the ratio of the up-quark and down-quark masses. Likewise, the coefficients CairN and 
Ca-n for such an axion are 

^'^'^^ = 2V2(l + z) ' ^^"^ ^ W^) ' ^^-^"^ 

where /jr ~ 93 MeV is the pion decay constant and ttItt w 135.0 MeV is the neutral pion mass. 

Before concluding this review, we note that the effective coupling coefficients c^, Cg, and Ci appearing in Eq. (|2.2p 
are highly mo del- dependent. They need not be 0(1), and in many theories any of them may vanish outright. Indeed 
it has been argued _9] that the existence of axions and axion-like fields which couple to electromagnetism but not to 
SU{3) color is a generic feature of certain extensions of the SM, including string theory. In assessing the constraints 
on our bulk-axion DDM model, we shall therefore focus primarily on a "photonic" axion of this sort — i.e., a general 
axion with Cg = and Cy ^ 0. However, we shall also discuss how such phenomenological constraints are modified 
in the case of a so-called "hadronic" axion with non- vanishing values for both Cg and Cj. We note that additional 
subtleties arise in this latter case, due to non-trivial mixings between the a„ and other pseudoscalars present in the 
theory which also necessarily couple to G't^G'^^'^ . These include hadrons such as 7r° and 77, as well as any other axions 
in the theory which play a role in addressing the strong-CP problem [J, |5|. In discussing constraints on hadronic 
axions, we shall implicitly assume that the full mass-squared matrix for the theory is such that the relationship 
between the a„ and the mass eigenstates a\ defined in Eq. ()2.14p is not significantly disturbed. Indeed, given the 
inherently large number of independent scales and couplings that emerge in scenarios involving multiple axions and 
other pseudoscalars, this is not an unreasonable assumption; moreover, it is straightforward to show in a general way 
that these favorable conditions can always be arranged for certain sets of axion and pseudo-scalar mixings. Such an 
assumption thereby enables us to perform our phenomenological analysis in a model-independent way. 

III. AXION PRODUCTION IN THE EARLY UNIVERSE 

Axions and axion-like fields can be produced via a number of different mechanisms in the early universe. For 
example, these particles can be produced thermally, via their interactions with the SM fields in the radiation bath. 
In addition, a number of non-thermal mechanisms exist through which a sizable population of axions also may be 
generated. These include production via vacuum misalignment, production from the decays of cosmic strings and 
other topological defects, and production from the out-of-equilibrium decays of other, heavier fields in the theory. 



This last mechanism is particularly relevant in the context of the DDM models, since, by assumption, the dark sector 
in such models involves large numbers of unstable fields with long lifetimes. Indeed, in the axion DDM model under 
consideration in this paper, a non-thermal population of any a\ may be produced via the decays of both heavier KK 
gravitons and other heavier oa- 

In Ref. [21, we focused on misalignment production as the primary mechanism responsible for establishing a cos- 
mological population of dark axions. In order for the results for the relic abundances il\ of the a\ obtained there 
to be valid, the contributions from all of the alternative production mechanisms mentioned above must be subdom- 
inant for each ax. Therefore, in this section, we examine each of the relevant axion-production mechanisms in turn, 
beginning with a brief review of the results for misalignment production itself. Since phenomenological constraints 
on scenarios involving large, flat extra dimensions prefer that the reheating temperature T^ associated with cosmic 
inflation be quite low [7|, we will hereafter operate within the context of a low-temperature-reheating (LTR) cos- 
mology with Tfj ~ 0{MeV). Within such a cosmological context and within the region of model-parameter space 
in which misalignment production yields a total relic abundance fltot comparable to the observed dark-matter relic 
abundance Ocdm, we demonstrate that the contributions to each Ux from all other production mechanisms are indeed 
subdominant. 

A. Axion Production from Vacuum Misalignment 

We begin our discussion of axion production in the early universe with a brief review of the misalignment mechanism 
and its implications for axion DDM models. (A more detailed discussion can be found in Ref. ^.) As we shall see, 
this mechanism turns out to be the dominant production mechanism for dark-matter axions in such models. 

At temperatures T ^ Aq, the only contributions to the axion mass-squared matrix are the contributions from 
the KK masses. Since these contributions to Ai^ are diagonal in the KK eigenbasis, no mixing occurs, and the KK 
eigenstates are the mass eigenstates of the theory. The potential for each a„ with n ^ is therefore non- vanishing, due 
to the presence of the KK masses, and is minimized at a„ = 0. However, the potential for the zero mode ag vanishes. 
In the absence of a potential for oq, there is no preferred vacuum expectation value (VEV) (oq) which minimizes 
V{ao). It therefore follows that immediately following the phase transition a.t T ^ fx , the universe comprises a set 
of domains, each with a different homogeneous background value for the axion field which may be expressed in terms 
of a "misalignment angle" 9 = {ao)/fx- This angle is generically expected to be 0(1) in any particular domain, 
but could also be smaller. We assume here that Hj < "iirfx, where Hi is the value of the Hubble parameter during 
inflation, and therefore that the value of 9 is uniform over our present Hubble volume. In this case, we find that 

{ao)=9fx, (a„)=0 for n^O. (3.1) 

Note that the above discussion is strictly valid only in the limit in which the Hubble volume is taken to infinity. In 
reality, the presence of a finite Hubble volume limits our ability to distinguish fields with wavelengths larger than the 
Hubble radius from true background values. Because of this ambiguity, all a„ for which n/R < Hj can also acquire 
0(1) background values after U{l)x breaking. In Sect. lIVLl we will analyze the phenomenological consequences of 
this effect in detail and derive conditions under which it can be safely neglected. As we shall demonstrate, it turns out 
that within our preferred region of parameter space, these conditions involve only mild restrictions on the cosmological 
context into which our model is embedded. We will therefore assume from this point forward that the (a„) in our 
model are given by Eq. (|3.1I) . 

At temperatures down to T ^ Ag, the (a„) remain fixed at these initial values. At lower temperatures, however, 
the situation changes as instanton effects generate a potential for the axion KK modes. Indeed, in the regime in which 
T <^ Ag and the brane mass engendered by this potential has attained the constant, low-temperature value mx given 
in Eq. (|2.12p , the time-evolution of each field ax is governed by an equation of the form 

dx+'^dx+Txax + X'^ax ^ 0, (3.2) 

where each dot denotes a time derivative, and where 

in radiation-dominated (RD) eras , , 

in matter-dominated (MD) eras . 

When A < 3H/2, the solution to this equation remains approximately constant. This implies that the energy density 
stored in ax scales approximately like vacuum energy during this epoch. However, at later times, when A > 3iJ/2, we 



see that a\ oscillates coherently around the minimum of its potential, with oscillations damped by a "friction" term 
with coefficient 3H + r>. During this latter epoch, the energy density stored in a\ scales like massive matter. 

At temperatures T ~ Aq, the evolution of a\ depends more sensitively on the explicit time-dependence of the brane 
mass mx(t). In what follows, we adopt a "rapid-turn-on" approximation, in which the instanton potential is assumed 
to turn on instantaneously at t — Iq, where to is the time at which the confining transition for the gauge group G 
occurs. In this approximation, mx(t) takes the form of a Heaviside step function: 



mx{t) = mxQ{t-tG) 



(3.4) 



In this approximation, the a„ remain fixed at the initial values given in Eq. (13.11) so long as t < Iq. At t = tQ, the 
brane mass immediately assumes its constant, late-time value mx- Since only oq is populated immediately prior to 
the phase transition at ic, each of the ax initially acquires a background value proportional to its overlap with oq: 



(axita)) = ejxAx , {axita)) = . 
Subsequently, after the ax have been populated, each begins oscillating at a characteristic time scale 

tx 



max{^,iG} , 



(3.5) 



(3.6) 



where kx is the value of k corresponding to the epoch during which this oscillation begins. At late times t ^ tx, when 
these oscillations become rapid compared to the rate of change of (ax) and the virial approximation is therefore valid, 
one finds that the energy density px stored in each mode is given by 



Pxit) = ieViA^^^A ( J 



-T^it-tG) 



(3.7) 



during the epoch in which the oscillation began. Computing px during subsequent epochs is simply a matter of 
applying Eq. (13.71) iteratively with the appropriate boundary conditions at each transition point. Consequently, in 
the LTR cosmology, we have [^ 
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where i^u denotes the reheating time — i.e., the time at which T — Trh, and the universe transitions from an initial 
epoch of matter domination by the coherent oscillations of the inflaton field to the usual radiation-dominated era. 

Given the energy-density expression in p.8p . it is straightforward to obtain the relic abundance fi^ = px/ Petit for 
each ax, where pcrit = iMpH^. For the heavier modes in the tower, for which tx = tQ, one finds 
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For the modes in the tower for which tx > tc, the corresponding result is 
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The total contribution fitot to the dark-matter rehc abundance from the axion tower is simply the sum over these 
individual contributions. While the generic behavior of fitot as a function of fx, Mc, and Aq is somewhat complicated, 
simple analytical results can be obtained in certain limiting cases of physical importance. For example, let us consider 
the limit in which t\ — Iq for all modes in the tower and Hj is sufficiently large that none of the a\ which would 
otherwise contribute significantly to iltot begin oscillating before the end of inflation. In this limit, all of the Q,\ take 
the form given in Eq. \?>M . and one finds that the present-day value of fitot, here denoted H^^t, is given by the simple 
closed-form expression [2] 

^i^-^,i^oif{^£)e^v£..{^y . (3.1) 

In the opposite limit, when all of modes which contribute significantly toward $7*^^ begin oscillating at t\ > Iq and 
have oscillation-onset times which depend on A and are therefore staggered in time, U,\ is given by Eq. p.lOp for all 
ax. In this hmit, the expression for ri^^^ reduces to [4I 

2 



---107) IT^/"- <"'' 



n;,„ ^l "-m f^^") 



The preferred region of parameter space from the perspective of dark-matter phenomenology is that within which 
rijot represents an 0(1) fraction of the dark-matter relic abundance inferred from WMAP data flO ]: 

^cdm/i^ = 0.1131 ± 0.0034 , (3.13) 

where h « 0.72 is the Hubble constant. From a dynamical dark-matter perspective, it is also preferable that the 
full axion tower contribute meaningfully to 57;*^^. For an 0(1) value of the misalignment angle and a reheating 
temperature within the preferred range T^ ~ 4 — 30 MeV for theories with large extra dimensions, one finds [2| that 
these two conditions are realized for fx ~ 10^^ — 10^^ GeV and Aq ~ 10^ — 10^ GeV, provided that Mc is small 
enough that y ^ 1- Within this region of parameter space, the t\ of all a\ which contribute meaningfully toward 
rij*Q). are staggered in time, and therefore the lighter modes yield a proportionally greater contribution to that total 
abundance. We will often focus our attention on this particular region of parameter space when discussing constraints 
on axion DDM models. 

B. Axion Production from Particle Decays 

Another mechanism by which a non-thermal population of relic particles may be generated in the early universe 
is through the decays of heavier, unstable relics. In scenarios involving extra dimensions, these relics include the 
higher KK modes of any fields which propagate within at least some subspace of the extra-dimensional bulk. For 
example, since the graviton field necessarily propagates throughout the entirety of the bulk, a population of unstable 
KK gravitons is a generic feature of all such scenarios. In the minimal bulk-axion DDM model under consideration 
here, the unstable relics whose decays can serve as a source for any given a\ include these KK gravitons as well as 
other, heavier a\. Moreover, since these fields span a broad range of masses from the sub-eV to multi-TeV scale 
and beyond, one would expect the population of axions produced by their collective decays to possess by a highly 
non-trivial phase-space distribution. However, as we shall demonstrate below, the total contribution F^^^-* to the 
decay rate of any a\ from intra- ensemble decays (i.e., decays to final states which include one or more dark-sector 
fields in addition to any visible-sector fields that might also be present) is far smaller than that from decays to final 
states involving visible-sector fields alone. That the total branching fraction for intra-cnscmble decays is negligible 
suggests that the population of axions produced by such decays will, in general, be quite small. Thus, provided the 
initial abundances of the a\ are set by some mechanism such as vacuum misalignment for which the 17 a of the heavier 
a\ are initially similar to or smaller than those of the light fields, it is reasonable to assume that the contributions 
from intra-ensemble decays are subleading and may therefore be safely neglected. 

One class of processes which contribute to F^^^ are those which arise due to the axion self-interactions implied 
by the final term in Eq. (|2.17p . The leading such contribution comes from three-body decay processes of the form 
o-x ^ a-\xa-\20-\3- An upper bound on the total contribution V{a\ — > 3a) to the decay width of a given a\ from all 
kincmatically allowed decays of this form was derived in Ref . [2il : 



It was also shown in Ref. [4I that the partial width of the a\ to a pair of photons is given by 

r(aA^77) = G^(PAa)2^, (3.15) 

JX 

with G~f = c^a^/2567r'', where a = e^/47r is the fine-structure constant. Within the preferred region of parameter 

space discussed above, in which fx ^ 10^^ — 10"^^ GeV and Kq > 10^ — 10^ GeV, we see that Fa (a — >■ 3a) is negligible 
compared to T(ax — > 77). It then follows that Tx{a — >■ 3a) represents a vanishingly small contribution to the total 
decay width Tx of any ax in any theory with an 0(1) value of c^. We therefore conclude that decays of the form 
o-x — >■ o-XiaxzO-Xs do not play a significant role in the phenomenology of realistic bulk-axion models of dynamical dark 
matter. 

In addition to these decays, however, an additional set of decay channels — those involving lighter graviton or 
radion fields in the final state — are also open to the ax- In order to assess whether such decay channels are capable 
of yielding a significant contribution to the relic abundance of any of the ax, we begin by identifying the relevant 
interactions among the modes in the KK graviton and axion towers. Since we are considering the case of a flat extra 
dimension and assuming fluctuations of the metric to be small, it is justifled to work in the regime of linearized gravity. 
The relevant term in the five-dimensional action is therefore 



d^x / dy——TMNh^'^ , (3.16) 

where Tmn is the stress-energy tensor, and Hmn is the metric perturbation defined according to the relation 

2 

gMN = VMN H TJT^hMN ■ (3-17) 

The piece of the stress-energy tensor which involves the five-dimensional axion field a includes both a bulk contri- 
bution and a contribution arising from terms in the interaction Lagrangian which involve interactions of the axion 
with the brane-localized fields of the SM. The bulk contribution is given by 

Tmn = dMadNa-^miNidpad^'a) . (3.18) 

Upon KK decomposition, this contribution, when coupled to Hmn as in Eq. p.l6l) , gives rise to three-point interactions 
between a KK graviton or radion field and a pair of ax- These interactions lead to decays of the form ax — > G}!^ ay, 
where G/^™ denotes a KK graviton with KK mode number m. In the absence of an instanton-induced brane mass 
term mx for the axion field (i.e., in the mx — > limit, in which all mixing between the axion KK modes vanishes 
and ax — > a„), KK-momentum conservation would imply that only a single, marginal decay channel would exist 
for each ax. Hence the contribution to Tx from such decays can be neglected. However, the instanton contribution 
to the axion mass-squared matrix violates KK-momentum conservation, and therefore, despite the fact that these 
axion-axion-graviton interactions are Planck-suppressed, they can still potentially contribute significantly to Tx, due 
to the large number of modes into which each ax can decay. 
By contrast, the brane-localized contribution, which is given by 



fx 



^ c r 

2 Yl ~TJ2 [i^^i°')^^l'^l^^i + ('9^'a)1/'i7M7Vi - '2Vt^i^{dpa)ij,jP-f^iJr 



C7e^ . V.f. r. o „ f^oar. \ . ^9^ 






327r2/^/^ \ ■■ ' ■■ ■ / 327r2/^ 



(3.19) 



leads to four-, five-, and six-point interactions between the graviton field, the ax, and the various SM fields. These 
interactions take the same form as those discussed in Section IIIIBI save that each vertex involves the coupling of 
an additional KK graviton and is suppressed, relative to the corresponding interaction involving the axion and SM 
fields alone, by an additional factor of Mp. The rates for such interactions will therefore always be much smaller 
than those calculated in Section IIIIBI Indeed, even the total contribution to the decay rate of a given ax from such 
processes, summed over graviton KK modes, will still be suppressed by a factor of roughly M5, where M5 denotes the 
five-dimensional Planck scale, relative to the contribution from decays to SM fields alone. It is therefore sufficient, 
at least for our present purposes, to neglect T|j™° and to focus solely on the interactions arising from the bulk 
contribution T^f^^. 



10 

We begin our analysis of axion-axion-graviton interactions by expanding the five-dimensional axion field, as well as 
the various components /i^^, /i^s, and /155 of the metric perturbation Hmn, in terms of KK modes. The mode expansion 
of the axion field for the orbifold compactification considered here was given in Eq. ()2.3|) : the mode expansions of /i^^, 
/ip5, and /155 are analogously given by 






00 



fmy\ 
,. cos (—J 



'^'^vhn^ 



^m/i ,T sin I 1 



}i5 



00 



\ R J 






m=0 



\ R J 



(3.20) 



Note in particular that ^1^5 must be odd with respect to the parity transformation ^5 — ^ —x^. Upon substituting 
these KK-mode decompositions into the linearized-gravity action given in Eq. p.l6p and integrating over y, we find 
that the terms in the effective, four-dimensional interaction Lagrangian which govern the interactions between the 
graviton and axion KK modes consist of the following three contributions: 



27tR 






M, 



3/2 



dy- E 



?Vn?ViTp , (m) 
AMp ^^ 



(2(9^a(")a''a(P) - Ty^'^a^a^'^^apO^P)) A 



+ 

mnp 



-f,-(^)a(")a(^)A_^ 



"55^bulkj„, _ Y^ 'm'n'p,(m) 



M, 



3/2 



■'-dy= E 



m,n,p— 

00 00 



AMp 55 



9''a(")apa(P)A+„^ + (^) a(")a(P)A-„^ 






M, 



m— 1 n,p— 



(3.21) 



where 



A± 

rnnp 



— p ~r OjYip — 'fi Zt 0772^7T,-^p -f" (^772__72_p 



(3.22) 



For the purposes of computing Feynman diagrams, it is convenient to work in the unitary gauge, in which the h 



(m) 
7^5 



and /15™ fields with m > are set to zero by the five-dimensional gauge transformations guN "> gMN+dAi^N + dN^M, 
where cm is the gauge parameter. In this gauge, the contributions in the second and third line of Eq. p. 211) vanish 
(save for the interactions between the axion KK modes and the radion field /155 , which will be discussed in due time), 
and the physical, gauge-invariant graviton fields 






R 






R^ 






(3.23) 



reduce to /i}i™ for all m > 0. The relevant part of the effective Lagrangian consequently reduces to 



C 



(m>0) 



00 CxD 



EV^ In 

^—^ O. /o 



ht^ 



— 1 n,p— ^ ^ 



(2a^a(")a''a^P) - ?7^'^a''a(")apa(P)) A+„p + if (^) a^^^a^P^A 



(3.24) 



The expression in Eq. p. 241) can be rewritten in terms of the mass eigenstates ax via the mixing matrix U\n in 
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Eq. dHH). The result is 



oo oo 



L 



(m>0) 



oo oo 



m=ln,p=OA,A' V P L 



A 



mnp 



E E E ^^'^^"^^"A (2a''aA9'^av - v'-'d'^a.d.a,,) 



7n=l n^O A, A 



(' 



X I ''n-mf^„_m X' 



^n+mt^„+m^A/ T" ' rn-n'-J rn-n,X 



rm-nUl_„y + r_n-rnUi„_,^yj + V^" ^a^Gx, 



(n - m){rn^,nUl_^^y + rm-„C/^_„^;^,) + (n + m)(r„+„C/^^„^^, + r_„_™[/l„_.,„^;,,) ^ , (3.25) 



where in going from the first equality to the second we have exploited the Kronecker deltas in A^„ to evaluate the 
sum over p. It should be noted that in the notation employed in the above expression, Uj^^^ = by definition for 
n < 0. The sum over n in Eq. (|3.25p can also be performed analytically, and the resulting, final expression for the 
Lagrangian in terms of the ax is found to be 



oo _. r 



[f:^,2V2Mp 



m=l A, A 



(3.26) 



where the coefficients C^^Xx' ^^"^ ^mAA' ^^^ given by 



a 



(1) 

lAA' 



m^y^yX'^AxAy 



m-^y^ - 2m2y2(A2 + A'2) + (A2 - W^f 

~'2 V2r^2,,2/T2 , V2^ _ /T2 _ V2^2^ 



a 



(2) _ A\^\'^[m^y^{X^ + X'^) - (A^ - X'^Y]AxA 



lAA' 



y2]^^AyA _ 2TO2y2(A2 + A'2) + (A2 - A'2)2] 



(3.27) 



From the interaction Lagrangian in Eq. p.26p . it is straightforward to obtain the Feynman rule for the graviton- 
axion-axion interaction vertex in the unitary gauge: 



\/2Mf 



{kl^,k2l, + fci^A:2M - Vt^vki ■ fcj) C)^\x' " ^m^^^^^'^aa' 



Using this vertex rule along with the graviton-polarization sum rule given in Refs. [Ill Il2l|. we find 

(^i^) ^2 



|A^(a,^/i(™)av)| 



12M2(mMc 



A" + {mMcf + A'4 - 2(mM,)^A^ - 2(mAf,)^A'^ - 2A^A 



for A' < A. Consequently, the partial width of oa from such a decay is 

r(aA ^ h^;^)ay) = -^^(x'AxfCx'^Ayf [a^ + [mAQ^ + A'^ - 2{mM,f{X^ + A'^) - 2A2A'^ 



1/2 



(3.28) 



(3.29) 



(m). 



Once again, in order to obtain the full contribution T{ax -^ h^^a) to Tx from decays of the form ax — >■ hliu ay, it 
is necessary to sum over all combinations of final-state graviton and axion modes which are kinematically accessible. 
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As before, we will approximate the mode sums over both m and A' as integrals. This yields the result 



Am^iX-'Axf '■^ 



SttX^McM^ jao 



(A-A')/A/c _ 

dX' I dm{X'^Ayf 





1/2 



A* + (mil/4)4 + A'4 - 2(mMe)2(A2 + A'^) - 2A2a'2) 

(X-Xr 






dX'{X'^Ayf{X + X') 



{X' + X")E 



{X + X'Y 



2XX'K 



l\2 



(A-AQ 
(A + A')= 



where K{x) and E{x) denote the complete elliptic integrals of the first and second kind, respectively: 



K{x) 



7r/2 

^/l 



dO 



x^ sin I 



E{x) = 



f deVi^ 

Jo 



x"^ sin I 



, (3.30) 



(3.31) 



In order to compare T{a\ — > h^^a) to the rate for ax decays to SM fields, we numerically integrate Eq. (I3.30p 
over A' and compare the resulting expression to the decay rate T{a\ — > 77) to photon pairs. In Fig. [1] we plot the 
ratio T{a\ — > hfii,a)/T{a\ — >■ 77) as a function of A for a variety of different choices of fx- In each case, we have 
set Kq = 1 TeV, Mc = 10~^^ GeV, and £, = go — ^- It is evident from this plot that only for values of fx above 
the preferred range fx ^ 10^^ — 10^*^ GeV does the decay rate for ax — > h^i^ay become similar in magnitude to the 
rate for axion decays into brane fields. Indeed, for values fx within this preferred range, T{a\ — > h^^a)/T{a\ -> 77) 
never exceeds 0.06, even for the lightest modes in the tower. Furthermore, for values of fx of this magnitude, the 
lifetimes for all a\ light enough to have T{ax — )■ h^,^a)/T{a\ -^ 77) near this maximal value are parametrically larger 
than the present age of the universe, even once the additional contribution to F^ from a\ — )■ hfj^^ay decays is taken 
into account. Consequently, the decays of such fields are not cosmologically relevant, and since the branching fraction 
for all other, heavier a\ into final states involving KK gravitons is utterly negligible. We therefore conclude that 
intra-ensemble decays are not phenomenologically relevant for bulk-axion models of dynamical dark matter. 

Up to this point, we have focused chiefly on the effect of the tensor KK modes of the higher-dimensional graviton 
field on axion production in the early universe. However, we have yet to address the effect of graviscalars such 
as the radion on axion production. Since our minimal DDM model involves only a single extra dimension, only a 
single physical graviscalar mode (proportional to /155 ) appears in the theory. Furthermore, while the masses of the 

hjl^ are dictated by the compactification geometry alone, the mass of this radion field depends on the details of the 
mechanism through which the radius of the extra dimension is stabilized, and is consequently highly model-dependent. 
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FIG. 1: The ratio r{ax - 
Ag = 1 TeV and C = SG 



hf_iua)/T{a\ — s- 77), shown as a function of A for several different values of fx- Here we have set 
■ 1, and we have taken the compactification scale to be Mc = 10~^^ GeV. It is clear from this plot 



that this ratio is safely below unity for fx within our preferred region 10 — 10 GeV. 
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In this paper, we assume that the physical radion field is sufBciently heavy so as not to play a significant role in the 
decay phenomenology of the light ax fields which contribute significantly to rij^^. Nevertheless, we note that in 
scenarios which involve multiple extra dimensions of comparable size, or scenarios in which a specific model for radius 
stabilization is invoked, graviscalars may play a more significant role in the phenomenology of the dark sector. 

C. Axion Production from Cosmic Strings 

A population of cold axions can also be generated by the decays of topological defects. In our axion DDM model, 
this includes decays of the cosmic strings associated with the breaking of the global U{l)x symmetry. Such decays are 
relevant in situations in which this symmetry remains unbroken until after inflation, i.e., Hj > 2T:fx- By contrast, 
in situations in which Hi < Stt/x and the U{l)x is spontaneously broken prior to the infiationary epoch, cosmic 
strings and other topological defects are washed out by the rapid expansion of the universe during cosmic infiation. 
Consequently, in this latter case, axion production from the decays of cosmic strings can safely be ignored. 

In this paper, we are primarily interested in high values of fx ^ 10^^ — 10^*^ GeV, as these values characterize our 
preferred region of parameter space. Likewise, we will primarily be interested in relatively low values of Hj, which 
may be realized naturally in the LTR cosmology. For this reason we shall assume that Hj < li^jx in what follows. 
We see, then, that no significant population of axions is produced by cosmic-string decay. 

D. Axion Production from the Thermal Bath 

Another mechanism through which a relic population of axions may be produced in the early universe is a thermal 
one: via their interactions with the SM fields in the radiation bath. Unlike the axion population generated by vacuum 
misalignment, which is characterized by a highly non-thermal velocity distribution (essentially that of a Bose-Einstein 
condensate) and is therefore by nature cold, this population is characterized by a thermal velocity distribution. Indeed, 
the properties of a thermal population of axions can differ substantially from that of a population of axions generated 
via misalignment production. 

A number of processes contribute to thermal axion production in the early universe, and the processes which are the 
most relevant for the production of standard axions dominate for each a\ in this scenario as well. Among hadronic 
processes, which play an important role in axion production when Cg is non-vanishing, 97 — ^ qa\ and qg — > qa\ 
dominate for T > Aqcd, while pion-axion conversion off nuclei (including all processes of the form Nix — > N'a\, 
where N,N' — {n,p} and it denotes either a charged or neutral pion) and the purely pionic process tttt — > ttoa 
dominate at lower temperatures. The rate for the high-temperature process is [l3| 



--^raA^A^ 



Tiqj^qax) - ^^;^— -^ In 



647r5/ 



"^ -f- 0.406 



(3.32) 



where rUg is the plasma mass for the gluon, given in terms of the effective number of quark flavors Nf at temperature 
Tby 

^giT) = ^y^3 + iV//2. (3.33) 

Likewise, the rates for pion-conversion off nuclei and pionic production are well estimated by the expressions |14l4l6| 



7^7/2 3/2,4 ,2 -mjv/T r 

TiNn^N'a,) = ^ /' 1.64(5Cl + 5^7^^ + 2C,„C,p) + GCLat / 

6C(3)(27r)5/-^/^/^ L J Jo 

3C(3)T5CLA^A^ r r dx.dx^xlxl f' [s - mlf{hs~2ml) 



dx,^^ 
e^i - 1 



1024^7/2, /2 J^ J^ yiy2{ey^-l){ev--l)j-i " s^T^ 

where, once again, C,{x) denotes the Riemann zeta function, and the effective coupling coefficients Cap, Can, Cair, and 
Ca-nN are given in Eqs. (|2.19p and (|2.20l) . Since these processes are mediated by strong interactions, they tend to 
dominate the production rate for a hadronic axion at temperatures T > 100 MeV, at which the number densities of 
pions and other hadronic species are unsuppressed. 

In addition to these hadronic processes, there are several process involving the interactions between the a\ and the 
e and photon fields which contribute to the axion production rate, and indeed dominate that rate at temperatures 
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T ^ Aqcd- The first of tfiese is tfie inverse-decay process 77 —^ a^, tlie rate for which is given by 

r(77-«.) = 2^!^ie^i.,fA) (3.35) 

where Ki{x) and K2{x) respectively denote the Bessel function of the first and second kind, and G-y = a^c^/2567r^. 
Another is the Primakoff process 67 — > e a. For T, ?7ie 3> A, the rate for this process is weh approximated by [17[ 



s 2 



rprim(e±7 -> e±aA) = TT^^'^a 

192C(3)/1 



rp2 

In ( ^r 1 + 0.8194 



(3.36) 



where the plasma mass m^ of the photon is given by m^ = eT/3. In the approximation of vanishing chemical potential, 
the number density of electrons (plus positrons) Ue takes the well-known form 



3C(3), 



y3 J. > 



Tor 



~2^ 



3/2 



(3.37) 



Finally, if Ce 7^ in Eq. (|2.5p and the axion couples directly to the electron field, there can be an additional contribution 
to the e^7 — > e^a rate from a process akin to Compton scattering, but with an axion replacing the photon in the 
final state. The rate for this process can be estimated as |18| : 

rco.p(eS -. e±a.) ^ '-^{X^A^f x J ^ "^ ^ '"^ (3.38) 
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In Fig. [21 we provide a pictorial comparison of the rates for the axion-production processes enumerated above 
as functions of temperature. The left panel shows the rates for the production of a single axion species a\ with 
A = 1 MeV in a scenario with fx = 10^^ GeV, M^ = 10^" GeV, and Ac = 1 TeV. As before, we have taken 
^ = 1 and set Cg = c^ = c,^ = 1. The red curve corresponds to the rate r(7r7r — >■ Tiax) for the pionic process; the 
orange curve to the rate rprim(e^7 — >■ e^a\) for the Primakoff process; the green curve to the rate r(77 — >■ a) for the 
inverse-decay process; the blue curve to the rate V{N'k -^ N'a\) for the pion-conversion process off nuclei; and the 
purple curve to the rate r(g7 — >■ qa\) for the quark-gluon process. As the hadron description of the theory is valid 
only for T < Aqcd, and likewise, the quark/giuon description is only valid for T > Aqcd, the rates r(7r7r — >■ ttoa), 
r(A^7r — > N'ax), and r{qj -^ qa^) are only defined on one side or the other of this scale. The yellow curve corresponds 
to the rate rcomp('''"''' ~^ '^o,\) for the Compton-like process for Cg = 1. For purposes of comparison, we also show 
the Hubble parameter as a function of T for two different cosmologies: the standard cosmology (black dashed curve), 
and an LTR cosmology with a reheating temperature Trh = 5 MeV (black dot-dashed curve). The value of A we 
have chosen here is well within the asymptotic regime for this choice of Ale and fx', hence the rates displayed here 
represent take essentially the maximal values possible for any a\ in the scenario. In the right panel of Fig. [2l we show, 
for the same choice of Mc and fx , the total contribution to the axion-production rate obtained by summing the rates 
for all ax for which X <T — i.e., those which will be kinematically accessible at a given temperature. 

The most salient lesson to draw from of Fig. [2] is that even after the contributions from all kinetically accessible ax 
states are included in the thermal axion-production rate, none of the relevant processes by which a thermal population 
of axions might be produced comes close to satisfying the T ^ H criterion. This implies that the ax, even when taken 
together, never attain thermal equilibrium with the plasma after inflation ends. Furthermore, these results also 
justify the claims made above, that the electron Primakoff process and inverse decays of the form 77 — >■ a are the 
most relevant axion-production processes for T < Aqcd, while hadronic processes dominate the axion-production 
rate for T > Aqcd- 

Let us now estimate the contribution to Utot from thermal axion production in the context of an LTR cosmology 
with a reheating temperature of Trh = 5 MeV. For concreteness, we focus on the case of a photonic axion with Cj = 1 
and Cg = Ci — for all i; however, the results for other coupling assignments should not differ drastically from those 
obtained here. We begin by noting that any contribution to Qx generated at temperatures T > Trh, *-e-, during the 
reheating phase, will be substantially diluted due to entropy production from inflaton decays. It is therefore legitimate 
to restrict our attention to axion production within the subsequent RD era. For a photonic axion, the processes which 
contribute to thermal axion production are inverse decays and e^7 — > e^a, the latter of which, since we are assuming 
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FIG. 2: A comparison of the rates associated with different axion-production processes in the early universe. Here we have 
taken Ma = 10"" GeV, /x = 10^^ GeV, Ac = 1 TeV, Trh = 5 MeV, and ^ = 1. The left panel shows the production rate 
for each process for an individual axion species a\ with A = 1 MeV (i.e., a value well within the asymptotic, large-A regime, 
where the rates are the least suppressed). The right panel shows the integrated production rate for each process, including 
contributions from all modes with \ < T . The most relevant processes for thermal axion production in this scenario are 



TTTT — > ttoa production (red curve), 67 — > e a\ via the Compton process (yellow curve) 



e^'^ 



e ax production via the 



Primakoff process (orange curve), inverse decays of the form 77 — >■ ax (green curve), production via the quark-gluon process 
qg — >■ qa (blue curve), and pion-production off nuclei (purple curve). It should be noted that the Compton process requires a 
non-zero electron-electron-axion coupling Ce, and that the curve shown here corresponds to the case in which Ce = 1. The value 
of the Hubble parameter as a function of T in both the Standard (black dashed curve) and LTR (black dash-dotted curve) 
cosmologies are also shown. 



Ce — 0, is dominated by the Primakoff process. The Bohzmann equation for the number density nx of each ax is 
therefore effectively given by 

nx + {3H + rx)nx - Cf (T) -K CP"'" (T) (3.39) 

for T < Trh, where C^"'^{T) and Cj^{T) are the contact terms associated with the electron-Primakoff and inverse- 
decay rates given in Eqs. p.36p and p. 351) . respectively. For T ^ A, me, these contact terms are well-approximated 
by the expressions 
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(3.40) 



where Ki{x) denotes the Bessel function of the first kind. To obtain a rough estimate of the relic abundance in 
situations in which either me or A is comparable to or greater than T, we modify the expression for Cx™^{T) given 
in Eq. (|3.40p by including an additional exponential factor e~('^+™=)/-^ to model the effect of Boltzmann suppression. 
From Eq. (j3.39p we estimate the relic abundance of axions produced by interactions with the SM particles in the 
thermal bath. To do so, we neglect the decay term and rewrite the resulting equation in terms of the quantity 
Y\ = nx/s, where s is the entropy density, in order to remove the Hubble term: 
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By numerically integrating this equation, we obtain an estimate of the thermal contribution 51 
rix of each ax at present time: 
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FIG. 3: Contributions to the individual mode abundances i}\ for a photonic axion from thermal production (dashed curves) 
and misahgnment production (sohd curves), plotted as functions of A for fx = 10^^ GeV (red curves), fx = lO'^^ GeV (orange 
curves), fx ~ W^"^ GeV (green curves), fx = 10^^ GeV (blue curves). The left panel displays the results for Ag = 1 TeV, 
while the right panel displays the results for Ag — 100 TeV. The other model parameters have been set to Mc = 10^^^ GeV, 
Trh = 5 MeV, ^ = gG^Cj = l. 



where g,,siT) is the number of interacting degrees of freedom present in the thermal bath at temperature T, and 
where k{T) is defined in Eq. (I3.3|) . The results of this integration are displayed in Fig. [3] In this figure, we compare 
the contributions to the relic abundance fix of a given ax from misalignment production and thermal production for 
a variety of different choices of the model parameters. 
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contribution 51 j^ from vacuum misalignment for reasonably heavy ax- Neither 51 
negligible compared with the 513^™"' contribution from the lighter modes. Indeed, summing over A to obtain the total 
thermal contribution 51^^^'^'^™'^ to the axion relic abundance at present time yields 3.8 x 10~^ < 5]^^^'^'^™ < 3.8 x 10~^. 

<C ^tot"' within the preferred region of parameter space for bulk-axion 
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We may therefore safely conclude that 51. 

models of dynamical dark matter, and that the population of ax generated by the misalignment mechanism dominates 

the relic density of the DDM ensemble. 

To summarize the results of this section, we have examined the primary mechanisms through which a cosmological 
population of DDM axions may be generated, including misalignment production, thermal production, and produc- 
tion by decaying relics. We have shown that within the preferred region of parameter space specified in Ref. [2], 

the contribution to the total present-day dark-matter relic abundance from misalignment production 51):""*' indeed 
dominates over the contributions from all other production mechanisms. This justifies the emphasis placed on mis- 
alignment production in Ref. [2] . Still, we note that although populations of axions produced via those other channels 
collectively represent a negligible fraction of 51tot, those populations can nevertheless play an important role in con- 
straining bulk-axion DDM models. For example, the thermal population of axions discussed above can still leave a 



significant imprint on the diffuse X-ray spectrum despite the small size of fl 



(the 



because 51 



(therm) 



« 51^'°") when 



A is large. We shall return to this point in Sect. II VG] where we will show that this imprint is nevertheless consistent 
with current observational limits. 



IV. PHENOMENOLOGICAL CONSTRAINTS ON DARK AXION ENSEMBLES 



In the previous section we characterized the various mechanisms which contribute to the generation of a cosmological 
population of relic axions in axion DDM models and compared the sizes of their contributions to 51^^^ . Given that 
this population constitutes the dark-matter ensemble in our axion DDM model, we now turn to examine the relevant 
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phenomenological, astrophysical, and cosmological constraints on that population of axions. As we shall see, some of 
these constraints pertain generically to any theory of dark matter, or to any theory containing late-decaying relics. 
Others are particular to models involving light, weakly-coupled fields. Still others pertain to theories with large extra 
dimensions in general, regardless of the presence or absence of a bulk axion field. 

As we have seen in Refs. [l|, y|> the properties of the dark- matter ensemble and its constituent fields in our bulk- 
axion DDM model are determined primarily by three parameters: the compactification scale Ale, the C/(l)x-breaking 
scale fx, and the confinement scale Aq for the gauge group G. Because these parameters play such a central role 
in characterizing the dark sector in our model, we shall seek to phrase our phenomenological constraints in terms of 
restrictions on Af^, fx, and A^ whenever possible. Of course, in addition to these primary parameters, a number 
of other ancillary quantities also have an impact on the phenomenology of our model, and thus are also constrained 
by data. These include the scales Hj and Trh associated with cosmic inflation, the coupling coefficients Cg, c-y, 
and Ci, and so forth. Generally speaking, these additional parameters play a subordinate role in determining the 
mass spectrum and relic abundances of the a\, and the values they take are typically far more mo del- dependent 
than Mc, fx, and Aq- Thus, while certain experimental and observational limits serve to constrain the values these 
additional parameters may take, it ultimately turns out to be possible to phrase the majority of constraints on our 
model as bounds on M^, fx, and A^. Indeed, as we shall see in Sect. |Vl most of the critical bounds can be expressed 
conveniently in this manner. We will also be interested in how these bounds constrain certain derived quantities of 
physical importance, such as the quantity y defined in Eq. (J2.12I) . which quantifies the amount of mixing that occurs 
across our DDM ensemble. 

A. Constraints from Background Geometry 

The first set of constraints we consider are those which apply generically to theories with extra dimensions, inde- 
pendently of the presence or properties of the bulk axion field whose KK excitations constitute the DDM ensemble 
in our model. These constraints arise primarily from experimental limits on the physical effects to which the tower 
of KK gravitons necessarily present in such theories gives rise. We will primarily focus here on scenarios involving n 
flat extra dimensions in which the fields of the SM are localized on a 3-brane, while gravity, as always, necessarily 
propagates throughout the entirety of the D — {4 + n)-dimensional bulk. 

Perhaps the most significant and direct bound on Mc in theories with extra dimensions arises due to modifications 
of Newton's law at short distances as a consequence of KK-graviton exchange. The lack of evidence for any such effect 
at modified-gravity experiments to date implies constraints on the sizes and shapes of those extra dimensions. For 
the case of a single large, flat extra dimension, the current limit on the compactification scale from such experiments 
is [111 

Mc > 3.9 X 10^12 GeV . (4.1) 

This lower limit on the compactification scale is robust in the sense that even if there exist additional compact 
dimensions with radii r^ <^ l/Mc, this bound is essentially unaffected. For this reason, Eq. (|4.1I) turns out to 
represent the most significant constraint on the parameter space of bulk-axion DDM models from considerations 
which derive solely from the presence of extra dimensions. 

There also exist additional constraints on the compactification geometry which arise due the relationship between 
this scale, the effective four-dimensional Planck scale Mp, and the fundamental scale of quantum gravity Mo. These 
constraints are generally more sensitive to the details of the compactification scenario. In general, the fundamental 
scale AId is related to Mp by [g] 

M^ = T4M^+" , (4.2) 

where Vn is the volume of the n-dimensional manifold on which the extra dimensions are compactified. For the simple 
case in which this manifold is a flat, rectangular n-torus, the volume Vn is simply the product of (27rri) for each cycle 
of the torus. Assuming all radii are equal to a common radius r, we then have 

r-1 > 2nMS"' i^-^j . (4.3) 

Bounds on the scale Mp) appearing in the literature are frequently predicated on these assumptions. However, we 
emphasize that in situations in which the r^ are not all equal, or in which the compactiflcation geometry differs from 
that of a flat, rectangular n-torus, those bounds can be considerably modified. 
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Under the assumption that the compactification geometry resembles that on which Eq. (|4.3[) is predicated, one may 
derive constraints on M^,, r, or combinations of the two. For example, one class of constraints which arise in theories 
with extra dimensions are those implied by the non-observation of effects related to thermal KK graviton production 
in astrophysical sources such as stars [2U| and supernovae [2l|, 122] ■ A brief synopsis of the most relevant bounds in 
this class is given in Ref. (20| , all of which depend crucially on the fundamental quantum-gravity scale Me, . The most 
stringent of these constraints currently derives from limits on photoproduction and stellar heating by gravitationally 
trapped KK gravitons in the halos of neutron stars. Indeed, for a theory involving n extra dimensions with equal radii, 
one finds that for n = 2, the bound is r"! > 5.8 x 10"'^ GeV, while for n = 3, one finds r'^ > 3.8 x IQ-^" GeV |20|]. 

Collider data also place limits on r and Md in theories with extra dimensions. Searches for evidence of KK- 
graviton production in the monojet (i.e., j + I$t) channel have been performed by the ATLAS [23, [SO] and CMS |31j 
collaborations. The most recent ATLAS analysis [30], conducted with 1 fb""'^ of integrated luminosity, constrains 
Md > {3.16,2.50,2.15} TeV at 95% C.L. for n = {2,3,4} flat extra dimensions with equal radii. The most recent 
CMS analysis [3l|, conducted at a comparable integrated luminosity, yields the slightly more stringent constraint 
Md ^ {4.03,3.21,2.80} TeV at 95% C.L. for the corresponding values of n. Limits from searches for KK-graviton 
effects in the diphoton [32] and dimuon [335 channels at 36 pb~ and 39 pb^ of integrated luminosity, respectively, 
have also been derived by the CMS collaboration, but these are currently less stringent than the constraints from the 
j + I^T channel. 

It is important to realize that the aforementioned bounds on Mjj as a function of the compactification geometry do 
not necessarily translate directly into analogous bounds on fx for a given fx- Unlike the graviton field, the bulk axion 
field in our DDM model need not necessarily propagate throughout the entirety of the extra-dimensional volume, but 
may in principle also be confined to a (4 -I- 7ia)-dimensional subspace of that volume, where Ha < n. When this is the 
case, fx is related to fx by the generahzation of Eq. ()2.7p : 



fx = K./i+"" . (4.4) 

Note that this relationship differs from that which exists between Alp and Md because ria < n. In this paper, 
as in Ref. [3], we focus on the case in which the axion field propagates in a single extra dimension of radius i?, 
irrespective of the size, shape, or number of extra dimensions which compose the totality of the bulk. Accordingly, 
we define Mc = 1/ R to be the compactification scale associated with this particular extra dimension, and we shall 
use this notation throughout. In this paper, we are not aiming to set Mjj at or even near the TeV scale, since we 
are not attempting to solve the hierarchy problem, but rather to address the dark-matter problem. We will therefore 
assume that the structure of any additional bulk dimensions is such that phcnomenological constraints on Md and 
the associated compactification geometry are satisfied. Note, however, that the Newton's-law bound in Eq. (|4.ip does 
apply to Mc, as it applies to the compactification scale associated with any individual extra dimension. 

Another class of constraints on scenarios involving large extra dimensions applies to ancillary variables which charac- 
terize the cosmological context in which our model is situated. For example, the prediction of the observed abundances 
of the light elements via big-bang nucleosynthesis (BBN) is one of the greatest successes of the standard cosmology. 
Consistency with these predictions requires that effects stemming from the presence of these extra dimensions not 
disrupt BBN. Successful nucleosynthesis requires that the expansion rate of the universe during the BBN epoch, as 
quantified by the Hubble parameter H{T), must not deviate from its usual, four-dimensional value by more than 
around 10% [3]. In other words, there exists some temperature T„ > Tbbn ~ 1 MeV (usually dubbed the "normalcy 
temperature" in the literature) below which the radii of all extra dimensions are effectively fixed and the bulk is 
effectively empty of energy density. A variety of different considerations constrain T* , most of which are related to 
the potentially observable effects of KK-graviton dynamics in the early universe: 

• Interactions between the SM fields on the brane and the bulk graviton field result in a transfer of energy from 
the brane to the bulk, and a consequent cooling of the radiation bath on the brane. Substantial energy loss 
via this "evaporative cooling" mechanism would result in a modification of the expansion rate of the universe. 
At temperatures T < Tbbn, such a modification would distort the light-element abundances away from those 
predicted by standard BBN. Thus, the strength of the interactions between SM particles and excitations of the 
graviton field is constrained. 



• 



If the collective energy density associated with the graviton KK modes is substantial, that energy density could 
cause the universe to become matter-dominated too early. In extreme cases, it could even overdose the universe. 



• Late decays of KK gravitons could result in distortions of the abundances of light elements away from the values 
predicted by BBN [7], which accord well with the observed values for these abundances. Such decays could also 
result in significant entropy production. 
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• The relationship between the Hubble parameter H and the total energy density p of the universe is modified 
at early times in higher-dimensional scenarios, even when that energy density is overwhelmingly dominated by 
brane-localized states [23l - l25l |. Such a modification could have a substantial eS'ect on BBN as well. 

The constraints on T* implied by these considerations have been reckoned by a number of authors [7|, |26[ , and while 
the precise values of the bounds so derived again depend on the number, size, and shape of the extra dimensions, the 
value of M^), etc., the most stringent (which tend to come from limits on the late decays of the excited KK modes) 
generally tend to restrict T* to within the rough range 4 MeV '^T^ <2Q MeV [3] . 

One possibility for achieving such conditions is to posit that T* be identified with the reheating temperature Trh 
associated with a period of cosmic inflation initiated by an inflaton fleld which is localized on the same 3-brane as 
the SM fields. During such an inflationary epoch, any contributions to the energy density of the universe from bulk 
states which existed prior to the inflationary epoch (save for those which, like the contributions to p\ from vacuum 
misalignment, scale like vacuum energy) are inflated away. Furthermore, if the inflaton fleld decays primarily to 
other brane-localized states, no substantial population of bulk states is regenerated during the subsequent reheating 
phase. Thus, by adopting a LTR cosmology with a reheating temperature 4 MeV < Trh ^ 20 MeV, we thereby 
ensure that the relevant constraints related to KK-graviton production in the early universe are satisfled. We also 
note that a reheating temperature of Trh ^ 4 MeV is sufflcient to ensure that the thermal populations of the SM 
fields (and, in particular, the three neutrino species) required in standard BBN are generated by the thermal bath 
after reheating [13, [H. 

In summary, while stringent constraints exist on theories with large extra dimensions, these constraints can be 
satisfied by adopting an LTR cosmology with 4 MeV < Trh ^ 20 MeV and a compactification manifold for which 
the astrophysical bounds listed above may consistently be satisfied for a given choice of Md and fx- Since 0,1^^ is 
generated via non-thermal means in our bulk-axion model, as discussed in Sect. IIIIl such a low value of Trh is not 
an impediment to obtaining a dark- matter relic abundance rij^j « ^cdm- In fact, as shown in Ref. [2|, adopting an 
LTR cosmology is actually an asset in terms of generating a dark-matter relic abundance of the correct magnitude. 
Likewise, since the relationship between fx and fx need not be identical to the relationship between Mp and M/j, 
constraints which concern the effects of KK gravitons can be satisfied without imposing equally severe restrictions on 
the parameters fx, Mc, and Ag which govern the properties of the dark-matter ensemble. Indeed, the only significant 
model-independent constraint on these parameters turns out to be the constraint quoted in Eq. (14.11) from tests of 
Newton's law at short distances. 

B. Axion Production with Subsequent Detection: Helioscopes and Light Shining Through Walls 

We now address the constraints which relate directly to the phenomenological, astrophysical, and cosmological 
implications associated with the KK tower of axion fields which constitute the DDM ensemble in our model. We 
begin by discussing the limits derived from a wide variety of experiments designed to detect axions and axion-like 
particles via their interactions with the photon field. (For extensive reviews of these experiments, see Refs. |34ll35|.) 
To date, none of these experiments have seen any conclusive evidence for such particles, and the null results of these 
experiments therefore imply constraints on the effective couplings between such axion-like particles and the photon 
field. 

In order to determine how the results of the experiments listed above serve to constrain the parameter space of 
our bulk-axion DDM model, it is useful to divide those experiments into several broad classes, based on the sort 
of physical process each probes. One important class of experiments comprises those in which axions are produced 
via their interactions with the fields of the SM and then subsequently detected via those same interactions. These 
include helioscope experiments such as CAST [3g| and "light-shining-through- walls" (LSW) experiments such as BEV, 
GammaeV, and ALPS. Searches for coherent conversion of solar axions to X-ray photons in germanium and sodium- 
iodide crystals via Bragg diffraction which have been performed at experiments such as DAMA [33] , TEXONO [Sg , 
SOLAX [33], and COSME |4^] also fall into this category. The characteristic which distinguishes experiments in this 
class from others is that these experiments are affected by decoherence phenomena. Indeed, it has been observed Q 
that in theories with bulk axions, such phenomena result in a substantial suppression of the rate for any process 
involving the production and subsequent decay of axion modes relative to naive expectations. 

Let us briefly review the origin of this suppression by focusing on the interaction between the photon field and the 
axion KK modes given in Eq. (|2.5p . (The results for the coupling between these modes and the other SM fields are 
completely analogous.) We begin by defining a state 
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rnCn , (4.5) 
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which represents the particular hnear combination of KK eigenstates a„ that couples to any physics on the brane, 
such as Ff^^F^'^ or any pair of SM fields. Here N ~ fx/^c denotes the number of modes in the sum. Written in 
terms of a', the relevant term in the interaction Lagrangian becomes 



L,, 9 ^^a'F,.F^^ . (4.6) 



In other words, a' couples to the SM fields with a strength proportional to VN/jfx ~ I/Zjc- Consequently, the 
cross-section for any physical process which involves axions production via interactions with the SM fields followed 
by subsequent detection via the same sorts of interactions will take the form 



Ar2 
a{t) (X ^ X Pit) , (4.7) 
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where P{t) — |(a'(i)|a'(to))P is the probability for a state a' created at time to to be in the same state a' at time t. 
It can be shown that when N is large, P{t) is given by 
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where p is the initial momentum of the axion. 

At very early times, when t w tg, the cosine factor in P{t) is approximately unity for all values of A and A'. At 
such times, all of the terms in the sum appearing in the second term on the right side of Eq. (14.81) add coherently. 
As a result, this term, combined together with the first term, yields a factor on the order of A^^. However, as the 
system evolves, the cosine terms will no longer sum coherently, and a random-walk behavior ensues, according to 
which the two terms combine to yield a factor of 0{N) rather than of 0{N'^). The time scale td associated with this 
decoherence — or, more precisely, the scale at which P{t) = 0.1P{to) — is found to be [3] 

where y is defined in Eq. (|2.12p . Since td is clearly quite small for any combination of p and fx values of experimental 
relevance, any method of detecting axions which relies on their production and subsequent detection will feel the effect 
of this decoherence. By contrast, detection methods which rely on axion production without subsequent detection 
(such as missing-energy signals at colliders, energy dissipation from supernovae, etc.) or which probe for evidence of 
a cosmic population of relic axions (such as microwave-cavity experiments) will be unaffected by this phenomenon. 

The consequences of axion decoherence for physical processes in the decoherence regime are readily apparent. In 
this regime, as discussed above, the term in brackets in P{t) scales like A^ rather than A^^; hence any cross-section 
which takes the form given in Eq. (|4.7|) will scale with A^ ^ fx/M^ according to 

^11 . ^ 

^(t>^^) « _ (4.10) 

Jx ■' X 



In other words, such cross-sections are suppressed by an additional factor of A^ relative to the naive expectation 
obtained by setting fx — > fx in Eq- (|4.7p . Thus, due to the decoherence effect, any experimental bound on the 
effective coupling Ga-y-y of a single four-dimensional axion to the photon field which takes the form G^^^ < (GJJi^i')^ 
translates to a bound G^^^ < {G^^^)^/\/N for five-dimensional axion, rather than to G^^^ < {G™^^)^/N. Given the 
parametrization for Ga-y-y given in Eq. ()2.5p . we can phrase any such constraint as a bound on fx'- 
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Using Eq. (12. 7p . we may rewrite this constraint in the form 
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The most stringent limit from the class of experiments categorized above (i.e., those for which the phenomenon of 
decoherence is relevant) is currently the Ga-y^ < 8.8 x 10^^^ GeV~^ bound obtained by CAST [36|]. The most stringent 
hmit from crystalline detectors is the Ga-y-y ^ 1-7 x 10^^ GeV^ bound from DAMA [33|, and limits on Ga-y-y from 
LSW experiments are typically roughly three orders of magnitude higher than the CAST limit. The corresponding 
bound on fx from Eq. ()4.12|) is 

fx > (2.92 X 10«) 4/^ (^^) '^' GeV . (4.13) 

Note that even for M^ at the experimental lower limit given in Eq. (|4.1[) . the constraint in Eq. (|4.13p is satisfied as 
long as fx > 5.58 x lO^" GeV. 

C. Microwave-Cavity Experiments and Direct Detection of Dark-Matter Axions 

Another class of experiments which place constraints on the couplings of axions and axion-like fields to SM particles 
consists of those which involve the direct detection of a cosmological population of axions. The most sensitive 
experiments in this class are those associated with dedicated microwave-cavity detectors such as ADMX [4l| and 
C ARRACK [42l|. Detectors of this sort are used to search for the resonant conversion of dark- matter axions with 
mass ma to photons with energies E^ ss rria in the presence of a strong magnetic field. As a result, the observation 
of a signal at such a detector depends crucially on whether the mass of the axion in question lies within the range 
of photon energies probed. The axion mass range currently covered by ADMX spans only from 1.9 x 10~^^ GeV to 



3.5 x 10~^^ GeV |4l|, and the projected future mass sensitivity extends only as high as lO^^'^ GeV. Likewise, the 
projected sensitivity for CARRACK extends only as high as 3.5 x 10"^"' GeV. 

As discussed in Ref. Qj , the region of parameter space which is the most interesting from a DDM perspective is that 
within which y ^ 1 and mixing among the light axion KK modes is substantial, for it is this region within which the 
full tower contributes meaningfully to £7*^^ . Within this region of parameter space, the lightest mode in the tower has 
a mass Aq ~ Mc/2. Taken in conjunction with the bound on Mc from modified-gravity experiments given in Eq. (14.11) . 
this result implies that Aq > 1.5 x 10^^^ GeV in highly-mixed bulk-axion scenarios. The projected ranges for both 
ADMX and CARRACK lie well below this threshold for Aq. We therefore conclude that no meaningful constraints 
on bulk-axion DDM models can be derived from the results of these experiments. 

D. Axion Production w^ithout Subsequent Detection: Stars and Supernovae 

We now turn to examine an additional class of constraints on bulk-axion DDM models: those related to astrophysical 
processes in which the ax are produced through their interactions with the SM field, but never directly detected. 
Among the constraints in this class are limits on axions, moduli, and other light scalars derived from the non- 
observation of their would-be effects on the lifetimes, energy-loss rates, etc., of various astrophysical sources such as 
stars and supernovae. These effects include the following: 

• Axions and other light fields whose interactions with the particles of the SM are extremely weak and whose 
mean free paths are consequently extremely long can dissipate energy from stars extremely efficiently. Such 
dissipation can accelerate stellar cooling and result in observable alterations in stellar life cycles, including the 
life cycle of our own sun. 

• Similarly, such light fields can carry away a substantial fraction of the energy liberated by supernovae. Limits 
may therefore be placed on the strengths of these interactions from the non-observation of such effects for 
supernova SN1987A. 

• A diffuse population of long-lived axions or KK gravitons initially produced by stars and supernovae could decay 
at late times, distorting light-element abundances and producing an observable X-ray or 7-ray signal in the keV 
— MeV range or higher. No evidence for such a signal has been seen by EGRET, FERMI, HEAO, Chandra, 
COMPTEL, etc. 

As is well known, these considerations lead to some of the most stringent constraints on standard, four-dimensional 
QCD axions. We now turn to examine how these limits constrain the parameter space of generalized bulk-axion 
models. 
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The primary distinction between processes in which the presence of the a\ is ascertained by direct detection and 
those in which it is only inferred from an energy deficit is that in the latter class of processes, the a\ appear as 
particles in the asymptotic final state. Thus, the contributions from the individual ax to the overall event rate for any 
such a process add not at the amplitude level, but at the cross-section level. The decoherence phenomena discussed in 
Sect. IIV Bl are therefore irrelevant for such processes, and the total cross-section af^° for the production of "missing 
energy" in the form of a\ fields by any given physical process is simply the sum of the individual production cross- 
sections (T^'° for each axion species. Since the effective coupling between each ax and any pair of SM fields includes 
a factor X^Ax/fx from mass mixing, as indicated in Eq. (|2.17p . each of these individual production cross-sections 
scales as 



prod 



1 



fx 



(4.14) 



When it occurs, axion production will have a characteristic energy scale E^h determined by the surrounding envi- 
ronment. This energy scale may be associated, for example, with the temperature of a star or supernova core, or with 
the center-of-mass energy -^/s of a collider. Provided that Ech ^ Mc (an assumption valid for all physical contexts 
of relevance in bounding the large-extra-dimension scenarios considered here), it follows that A ^ Ech for a large 
number of ax- Such ax can be considered to be effectively massless as far as production kinematics is concerned, 
implying that to a very good approximation, a^" depends on A exclusively through the coupling-modification factor 
appearing in Eq. (|4.14p . (For those modes for which threshold effects are important, such an approximation will 
overestimate (T^'° and result in an overly conservative bound.) By contrast, a^^° will be effectively zero for those 
ax with masses A 3> -Ech in any thermal environment due to Boltzmann suppression, and will vanish outright in a 
non-thermal one. Therefore, it is reasonable to evaluate af^^ by taking any additional factors in Eq. (|4.14p to be 
essentially independent of A and by truncating the sum over modes at A ^ Ech- Thus, we find that 
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where the "effective" coupling H(_Ech) is given by 
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(4.16) 



Since the number of modes contributing to af^° is large by assumption, and since their masses are closely spaced, it 
is generally legitimate to approximate H(_Ech) by an integral 
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(4.17) 



The quantity H(i?ch) clearly plays a crucial role in the phenomenology of bulk-axion scenarios. It is therefore worth 
pausing a moment to examine in detail how H{Ech) depends on the physical scales fx, Mc, and Aq. A straightforward 
calculation shows that H(i5ch) has the parametric scaling behaviors 
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(4.18) 



The first case in Eq. (|4.18p corresponds to y <C 1, signaling a highly- mixed axion KK tower for which X'^Ax ^ 
X. By contrast, the second case corresponds to y 3> 1, signaling a relatively unmixed axion KK tower for which 
X^Ax ~ constant. These results for H(i?ch) are illustrated in the left panel of Fig. 0] for Ech — 30 MeV, a value which 
is physically meaningful in that it corresponds roughly to the core temperature of supernova SN1987A. Remarkably, 
we observe that H(i<^ch) experiences a suppression for y <^ 1. In other words, mixing within the axion KK tower acts to 
suppress the magnitude of the total production cross-section for processes in which the ax appear as missing energy. 
This is an important result, for it indicates that constraints on the parameter space of our bulk-axion DDM model 
derived from limits on axion production in stars, supernovae, colliders, etc., will be considerably weaker than one 
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FIG. 4: The dimensionless "effective coupling" factor b!(i?ch) defined in Eq. (|4.f6|l . shown as a function of the relevant scales 
Mc, fx, and Kg- In the left panel, we display curves of K(i5ch), each corresponding to a particular value of Ag and normalized 
to the value (-Ech/Mc)^'^ taken by K(ii'ch) in the absence of mixing, as a function of jx with fixed Mc = lO"'^^ GeV. It is 
readily apparent that the net effect of mixing within the KK axion tower is to significantly suppress this effective coupling, 
thereby loosening the corresponding production-cross-section constraints. In the right panel, we display curves showing the 
overall cross-section-suppression factor N^(-Bch)//x as a function of /x for fixed A = 1 TeV, each corresponding to a particular 
value of Afc- For each set of curves, we have taken ^ = po = 1, and have chosen _Ech = 30 MeV, which corresponds roughly to 
the core temperature Tsn of SN1987A. 



might expect from naive dimensional analysis. Moreover, this result applies more generally to any theory involving 
KK towers of scalar fields whose squared-mass matrix contains both brane-mass and KK-mass terms. 

In order to illustrate more explicitly the physical consequences of H(£'ch) in our bulk-axion DDM model, we likewise 
display the behavior of the overall scaling factor H^(i?ch)//x for erf™ in the right panel of Fig. |4l The results shown 
in this panel further illustrate a significant general property of this scaling factor: namely, that the cross-section is 
actually suppressed not only for large /x, but also for small /x, due to the parametric behavior of H(£^ch) described 
in Eq. (|4.18p . Thus, for any given choice of Aq and M^.^ there exists a maximum possible value for a\]^ , which is 
only attained at some particular value of jx- These results again illustrate the dramatic effect that H(_Ech) can have 



m suppressing cr. 



prod 



in our bulk-axion DDM model. 



Within the class of constraints from processes in which axions are produced but not subsequently detected, use of 
H(£'cii) allows us to translate experimental bounds on four-dimensional axion models into bounds on theories including 
towers of bulk scalars. The leading such bound is obtained from energy-loss limits from SN1987A. For a standard 
four-dimensional QCD axion, this bound is roughly [43| 



/a > 4 X 10^ GeV 



(4.19) 



By contrast, in the bulk-axion scenario under consideration here, each a\ light enough to be produced within the 
thermal environment of SN1987A can contribute to the overall energy-dissipation rate. Since the temperature Tsn 
associated with the supernova core is roughly 30 MeV, the appropriate modification of Eq. (|4.19p for a general axion 
which couples to hadrons with a coupling coefficient comparable in magnitude to that of a QCD axion is 



/x 



> 



(4 X 10* GeV) H(Tsn) 



(4.20) 



It then follows that in highly-mixed scenarios, this constraint can be significantly weaker than the corresponding 
constraint on KK-graviton production derived in Ref. 7], due to suppression by H(Tsn)- Indeed, the corresponding 
constraint on KK-graviton production is directly obtained by replacing ]x -^ ^Ip and H(Tsn) —^ ^sn/A/c in Eq. (|4.20|) . 
While the SN1987A bound is indeed one of the most stringent constraints on the QCD axion, it is not necessarily 
applicable for all general axions. This is because the bound quoted in Eq. (|4.19p is predicated on the assumption that 
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nucleon bremsstrahlung {N + N^^N + N + a) and other hadronic processes dominate the rate for the production 
of the hght scalar in question in the supernova core. This presupposes that the hght scalar couples to nuclei with a 
strength comparable to that of a QCD axion. If this is not the case, however, the constraints obtained from SN1987A 
energy-loss limits can differ considerably from the standard QCD-axion bound. For example, if from among the SM 
particles, the general axion couples only to the photon field, the dominant production processes will be e~7 — )■ ea\, 
p+7 — )• p^a\, and p^n — > p'^wyax. In this case, the considerably weaker bound 44] 

fx > (2.32 X 10^ GeV) c^ (4.21) 

is obtained for a four-dimensional field. Translating this result to the case of a KK tower of axions, as above, we find 
that 

fx > (2.32X 10^ GeV)c^H(TsN) . (4.22) 

Furthermore, in general axion models, c-y may not necessarily be of 0(1). In other words, the SN1987A constraint is 
sensitive to the U{l)x and SU{2) x U{1)y charges of the fields in the model, and is thus highly mo del- dependent. 

An analogous limit on fx can be derived from observations of the lifetimes of globular-cluster (GC) stars. The 
ambient temperatures Tqq of such objects are only 0(10 keV), so axion production primarily proceeds through the 
Primakoff processes 7-f e^ — )• a + e^ and j + nz — )• a + nz, where nz denotes a nucleus with atomic number Z. (Note 
that the dominant processes in this environment differ from the axion-nucleon- nucleon bremsstrahlung processes which 
dominate the axion-production rate in supernovae.) Such a bound will therefore arise for any general axion for which 
Cj ^ 1, regardless of whether or not it couples to the gluon field. The observation limit on axion production in GC 
stars is commonly phrased as an upper bound on the effective coupling Gay-y between a standard, four-dimensional 
axion (or any other similar particle) and a pair of photons, and the current bound is Ga^-y ^ 1 x 10~^° GeV~^ [45| . 
Since Tqc ~ 10 keV, the corresponding bound on fx is 

fx > (1.16 X 10" GeV) c^H(Tgc) • (4.23) 

Note that this constraint is independent of the SN1987A bounds, as it differs from the latter in two significant 
ways. First, because the relevant production process involves the coupling of the axion modes to photons rather than 
to nuclei, it depends on c-y alone and not on Cg. Second, since Tqc *^ 2sn, far fewer of the a a will be produced with 
any significant frequency in GC stars. Consequently, the enhancement factor from the sum over kinetically-accessible 
axion modes for GC stars is far smaller. 

Finally, bounds similar to those from SN1987A and GC stars can also be derived from the non-observation of effects 
related to axion production in other astrophysical sources, such as our own sun |46j . However, these bounds are found 
to be subleading in comparison with the SN1987A and GC-star constraints, essentially because they take place in far 
cooler environments, where the number of kinematically accessible modes is even further suppressed by the cutoff at 
Ech inherent in H{Ech)- 

E. Axion Production at Colliders 

We now consider the collider constraints applicable to our bulk-axion DDM model. Due to the highly suppressed 
couplings between the axion and the SM fields in standard four-dimensional axion models, collider data have virtually 
no relevance in constraining the parameter space of such models. Nevertheless, because of the huge multiplicity of 
light modes that arises in theories with light bulk fields in large extra dimensions, the net contribution to the event 
rates for certain processes from all of these modes taken together can potentially yield observable signals. For example, 
modes which are stable on collider time scales all appear as missing energy, and can lead to signals in channels such 
as pp ^ j + l^T and pp -^ j + Ifx- In addition, the heavier, more unstable modes which decay before exiting the 
detector can potentially give rise to additional signature patterns which may include displaced vertices. Indeed, we 
have already discussed in Sect. IIV Al how current limits from LHC data constrain the parameter space of one such 
bulk field — the higher-dimensional graviton — for which the monojet and monophoton channels mentioned above 
are of particular importance. Since the a\ in our bulk-axion model couple to the fields of the SM in much the same 
manner as KK gravitons, it is no surprise that the collider phenomenology of the a\ turns out to be quite similar to 
that of KK gravitons. 

We begin by discussing those signals which arise due to the combined effect of the a\ which are sufficiently long-lived 
so as to manifest themselves in a collider detector as missing energy. Collider processes in which the a\ appear as 
I^T are yet further examples of the class of processes discussed in the previous section in which axions are produced 
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but not subsequently detected. The net cross-section for any such process is therefore hkewise suppressed by axion 
mixing in the manner described in Eq. (|4.15p , with Ech given by the center-of-mass energy y/s of the coUider. 

Which specific channels are relevant for the discovery of a bulk axion at hadron colliders depends crucially on how 
the five-dimensional axion couples to the SM fields, and in particular on whether or not it couples appreciably to either 



light quarks or giuons. For a field with an 0{1) value of cither Cg 



or c, 



(where q — {u, d, s, c}), the principal discovery 



channel at both the Tevatron and the LHC is pp ^^ j + I$t, a- channel which is also one of the principal discovery 
channels for KK gravitons. Thus, in order to obtain a rough estimate of the constraints on the parameter space of 
our bulk-axion model from the null results of monojet searches, we translate the bound on the fundamental scale Md 
established by such searches into a bound on fx- The cross-section for KK-graviton production in association with 
a single jet at a hadron collider in a theory with n large, flat extra dimensions of equal length compactified on an 
n-torus, including contributions from all kinematically accessible modes, is roughly proportional to 



a 



3d(pp^j + G) ex 



27r 



1 



Ml^ 



(4.24) 



This implies that a bound of the form Mu > M™' 
our bulk-axion model of the form 



can be translated into a rough bound on the parameter space of 
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< 



X 



2tt 



(^/mm)„+2 



(4.25) 



where H(£'ch) is defined in Eq. ()4.16p . While this approximate bound does not take into account the differences in 
coupling structure between KK graviton and axion fields or the sum over polarizations for a massive graviton, it is 
sufhcient to obtain parametric estimates of the resulting constraints on our three fundamental parameters fx, Mc, 
and Ag- 

In Fig. m we indicate the rough bounds on the parameter space of our bulk-axion DDM model which can be derived 
in this manner, given a chosen value of Mf^™. The contours shown in this figure correspond to constraints of the form 
Md > A/]^'" for the illustrative values A/™'" = {1, 10, 100} TeV. 
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FIG. 5: Excluded regions of the {fx, Mc) parameter space of our DDM model in which the collider constraint in Eq. (I4.25|l is 
violated for Af]^'" = 1 TeV (green); M]j'" = 10 TeV (green and blue); and M]^'° = 100 TeV (green, blue, and purple). As Ac 
increases, we see that satisfying the collider constraints becomes increasingly easy, particularly for small fx- In each case, we 
have taken (, = Qc ~ I and assumed that the axion couples to at least one light, strongly-interacting SM particle with an 0(1) 
coupling coefficient Cg or Cq. 



We now compare these results to actual constraints on Mo from current experimental data and examine the 
projected LHC reach for our bulk-axion DDM model. The most stringent constraints from LHC data (which indeed 
come from the pp —>■ j + Ifr channel) were given in Sect. IIV Al Estimates of the future LHC reach for a theory with 
a single extra dimension are Mf^™ sa {14, 17} TeV at integrated luminosities dnt — {10, 100} fb~^, respectively j47| . 
Likewise, Tevatron data imply a limit Mf^™ k, 2.4 TeV for a theory with a single extra dimension 47[. Comparing 
these results to those in Fig. [SJ we see that current collider constraints, while quite stringent, do not significantly 
impact the preferred region of parameter space for our bulk-axion DDM model, even in cases in which the axion 
couples to one or more strongly-interacting SM fields with an 0(1) coupling coefficient. In such cases, since the most 
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stringent current LHC limits imply a bound of roughly M™'" « 1 TeV, the region of the parameter space of our model 
excluded by these limits roughly corresponds to the green shaded regions shown in Fig. [5] Since the green exclusion 
regions in this figure embody the most stringent such limits applicable to our DDM model, we shall take these to 
represent our collider constraints throughout the rest of this paper. However, we note that for photonic axions and 
other axion species which do not couple directly to quarks or gluons, the corresponding collider constraints (which 
arise from channels such as pp — > 7 + Ifr) are somewhat weaker. 

Before concluding this section, there is one important point which deserves emphasis. The collider processes we have 
been discussing thus far are those whose event rates receive their contributions from the low-lying modes in the tower 
— i.e., those a\ with lifetimes t\ > 10^^^ s. By contrast, those heavy ax with lifetimes t\ < 10~^^ s tend to decay 
to pairs of SM fields within the detector volume. The decays of such states can in principle give rise to an entirely 
different set of signature patterns. For example, a promptly decaying ax which couples to light quarks or gluons as 
well as photons would in principle contribute to event rates in the pp — >■ jjj and pp — > 77 + j channels. However, 
since the total event rate in these channels receives contributions from a broad spectrum of ax with different A, many 
event-selection techniques which are particularly useful in standard searches for new physics in these channels cannot 
be applied to a tower of decaying bulk axions. For example, since the set of decaying axions cannot characterized 
by a single, well-defined mass or cross-section, no identifiable peak can be expected to appear in the invariant-mass 
distribution for the decay products of the heavy axions. Such considerations render the results of standard searches 
for new physics in these channels inapplicable to our bulk-axion model — and indeed to DDM models in general. 
Moreover, they also likely render the identification of a conclusive signal of non-standard dark-matter physics in these 
channels particularly challenging. Nevertheless, the information that could potentially be revealed about the nature 
of the dark sector via such an identification is of sufficient magnitude and importance that an analysis of the discovery 
potential in these channels is an interesting topic for future study. 

F. Axion Decays and Distortions of the Cosmic Microwave Background Spectrum 

Up to this point, we have considered those phenomenological constraints on our DDM model which are related to 
the production of particles which compose our bulk-axion ensemble, both with and without their subsequent detection. 
By contrast, we now turn to discuss an entirely different set of phenomenological constraints, namely those which 
arise due to the potential decays of a pre-existing cosmological population of such particles. Indeed, such constraints 
emerge generically in all dark-matter scenarios in which the dark sector contains unstable, long-lived particles, and 
can be derived from observational limits on the physical consequences of the late decays of those particles. 

There are many considerations which can be used to place such limits on scenarios involving decaying dark-matter 
particles. For example, photons produced via the decays of such particles can yield observable distortions in the CMB 
spectrum; contribute to the diffuse extragalactic X-ray and gamma-ray backgrounds; upset BBN predictions for the 
primordial abundances of light elements; and result in unacceptably large entropy production during critical epochs in 
the history of the universe. Constraints on dark-matter candidates from considerations of this sort depend not only on 
the decay rate of the particle species in question, but also on the relic abundance of that species. For this reason, the 
constraints applicable to single-particle models of dark-matter are generally not directly applicable to models within 
the DDM framework. It is therefore necessary to revisit the observational limits on dark-matter decays within the 
context of our bulk-axion model of dynamical dark matter and assess how these limits constrain the parameter space 
of this model. 

In this section, we begin our analysis of the constraints on the late decays of the ax in our bulk-axion DDM 
model by examining observational limits on the distortions of the CMB which such decays can induce. The type of 
CMB distortion to which a late-decaying particle contributes depends on the time at which that particle decays. In 
the very early universe, photons produced by particle decays are brought into thermal and kinetic equilibrium with 
CMB photons via a number of processes. The dominant processes by which newly-produced photons can equilibrate 
thermally with CMB photons are double-Compton scattering (e~7 — >■ e~jj) and bremsstrahlung {e~X —>■ e~X 7, 
where X is an ion). However, once these processes freeze out, photons produced from ax decays are unable to 
thermally equilibrate with the radiation bath, resulting in the generation of a non-zero value for the pseudo-degeneracy 
parameter ^. The interaction rates for these processes are given by [48| 

Fbr « 1.57 X 10-36 (^1 _ ^^ {n^h'f/' (^^) '^' (^^) ''^' GeV , (4.26) 

where T^ow ~ 2.725 K is the present-day CMB temperature, Yp w 0.23 is the helium mass fraction, J7b ~ 0.044 is 



27 

the baryon density of the universe, and h « 0.72 is the Hubble constant. (Note that since z is quite large during 
the entirety of the relevant time frame, we have here approximated 1 + z « z.) Once these processes freeze out, in 
the sense that the rates given in Eq. (j4.26|) drop below the expansion rate H of the universe, photons produced by 
ax decay will no longer be able to attain thermal equilibrium with the CMB photons. Even after double-Compton 
scattering and bremsstrahlung effectively shut off, a number of photon-number-conserving interactions still serve to 
bring photons produced at even later times into kinetic equilibrium with the radiation bath. Dominant among these 
processes is elastic Compton scattering (e~7 — )■ e~j), which efficiently serves to bring photons produced by ax decays 
into kinetic equilibrium until a much later time isc ~ 9 x 10^ s, at which point this process too effectively freezes out. 
However, since elastic Compton scattering conserves photon number, it cannot similarly suffice to bring those photons 
into thermal equilibrium. As a result, CMB distortions in the form of a non-zero value for the pseudo-degeneracy 
parameter /i can be generated by ax decays during this epoch. In addition, after elastic Compton scattering freezes 
out, photons produced by ax decay achieve neither kinetic nor thermal equilibrium with the radiation bath. As a 
result, these photons no longer contribute the generation of /x, but instead contribute to the generation of a Compton 
y parameter (here denoted yc, so as to distinguish it from the ratio y = Mc/mx)- Finally, at t ^ 10^^ s, matter 
and radiation decouple, and any ax decays occurring after this point not affect the CMB, but instead simply persist 
as a contribution to the diffuse photon background. This last sort of contribution will be dealt with separately, in 
Sect. UVGl 

We thus see that axion decays have the potential to generate both a non-zero ^ and a non-zero yc- We can therefore 
establish constraints on our bulk-axion DDM model by calculating the theoretical predictions for these quantities in 
our model and comparing these predictions to observational data. 

We begin our analysis of CMB distortions from ax decays by addressing those decays which result in the generation 
of the pseudo-degeneracy parameter /x. In general, provided that the additional contribution 5p~j to the photon energy 
density p^ from the decay of the ax fields is small compared to the total p^, the time-evolution of p, can be described 
by the equation [H, ^ 






(4.27) 



Here Fdc smd Fbr are the interaction rates for double-Compton scattering and bremsstrahlung, respectively, and 
djia/dt denotes the differential contribution to /i from axion decay. For an arbitrary dfXa/dt, the solution to this 
differential equation takes the form 
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(4.28) 



where te » 1.69 x 10^ s is the time scale associated with electron-positron annihilation in the early universe, and where 
the quantities Cue a-nd Cbr are constants related to the double-Compton-scattering and bremsstrahlung rates Fdc 
and Fbr in Eq. (I4.26P by Fdc = Cuct^^^^ and Fbr = Cbr^"^^/^. Moreover, the differential contribution dp,a/dt to 
/i from axion decays is given by the standard expression for contributions due to the late injection of photons from a 
generic source: 
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(4.29) 



In general, the rate of change in the photon energy density is given by the Boltzmann equation for the evolution of 
p-y. In our bulk-axion DDM model, this equation includes a source term from each decaying state in the dark-matter 
ensemble. Thus, at late times, after all of the ax have already begun oscillating coherently and the contribution to 
Pj from inflaton decays can safely be neglected, we find that 
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AHp^ 



EBRf'^r.p. 



(4.30) 



where BRI '^ is the branching fraction of ax into a pair of photons. Note that the source term in the Boltzmann 



equation for p^ is simply a sum of the contributions from the various ax fields, 
relations 



Using Eq. (I4.30p . along with the 
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we can rewrite Eq. (J4.27I) in the form 
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(4.32) 



For the purpose of estabhshing a conservative bound, we focus here on the case of a purely photonic axion. As we 
saw in Sect. IIIIBl the contribution to Tx from intra-ensemble decays is neghgible for any a\ which decays on time 
scales relevant for the generation of CMB distortions. It is therefore justifiable to approximate Tx by the expression 
for r(a -^ 77) given in Eq. p.lSp and thus to take BRI^ '^ sa 1. Since the energy density px associated with each ax 
is given in Eq. p.7|) , we find that in this approximation, the first source term on the right side of Eq. (|4.32|) takes the 
form 
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where we have defined G-y is defined below Eq. p.lSp . The second term takes the same form, but with one factor of 
A fewer in the summand. 

In principle, one could evaluate this sum numerically at each moment in time, and then use these results to 
numerically solve Eq. (|4.32l) . However, we find that by making a few additional well- motivated approximations, we 
can obtain a closed-form, analytical result for dp,a/dt. We begin by dividing the tower into sections, based on the two 
criteria which determine the dependence of Vx and px on A. The first of these is whether the oscillation-onset time 
for a given ax is within the staggered regime (i.e., tx > tc), or the simultaneous turn-on regime (i.e., ^a = to). In the 
former case, tx depends on A according to Eq. l\'S.Q\i : in the latter case, tx is independent of A. The second pertinent 
criterion concerns the relationship between A and the quantity 
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(4.34) 



introduced in Ref. [2]. This quantity corresponds roughly to the transition point between the small- A regime, in 
which the ax are highly mixed, the large-A regime, in which mixing is negligible. Indeed, for A <C Atrans, we find 
that X^Ax ~ \/2 A/(l + 7r^/y^)^/^, while for A >• Atrans, we find that X^Ax ~ \/2. Given these criteria, our first 
approximation will be to replace X^Ax with its asymptotic large-A form for all A > Atrans, and with its asymptotic 
small-A form for all A < Atrans- Our second will be to approximate the sum over A by a set of source-term integrals 
Ii{m, n, a, /3, Amin, Amax), each corresponding to a different regime in the tower of modes characterized by a particular 
dependence of the integrand on A. These source-term integrals may be evaluated analytically by making use of the 
identity 
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(4.35) 



which is vahd for n > and any real values of ttz, a, and fS. Here r(5, x) denotes the incomplete gamma function: 

/•OO 

r{s,x) = / f-^e-^dt . (4.36) 



Employing the approximations discussed above, we find that the first source term on the right side of Eq. (j4.32p 
reduces to 
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TABLE I: Values of rrn, rn, Ui, and Pi which correspond to different regimes, labeled by the index i, in a generic axion tower, 
for use in Eqs. (|4.38|l and (|4.49p . The symbol kq denotes the specific value of k, as defined in Eq. p.3|) . which corresponds to 
ta- 



Inserting this result (and the analogous result for the second source term) into Eq. (|4.29p and using the fact that 
/i-type distortions are generated by decays occurring within the RD era, we obtain the result 
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(4.38) 



where the expressions for ai, /3i, tti^, and ni valid in each a\ regime are listed in Table HI Note that in obtaining this 
expression, we have assumed that the additional contributions to n-y and p-y due to the injection of photons from a\ 
are sufficiently small that these quantities can be approximated by the equilibrium expressions n'^ = 2C(3)T'^/7r^ and 
p°^ = 7r^r^/15. Furthermore, we have used the fact that the time frame during which CMB distortions to /i can 
arise lies entirely within the RD era. Obtaining a final result for the magnitude of /i-type distortions to the CMB 
engendered by the presence of a tower of decaying DDM axions is then simply a matter of substituting the result for 
dfia/dt in Eq. ()4.38p into Eq. (|4.28p and numerically evaluating the integral for a given choice of input parameters. 

The contribution to yc from the late decays of the a\ may be evaluated in much the same way as the corresponding 
contribution to /i. The decays which contribute to yc are those which occur during the window 9 x 10^ ^ ^ ^ 
1.2 X lO^'^ s, during which the rate Tec ~ H associated with elastic Compton scattering can no longer bring the 
photons from ax decay into kinetic equilibrium even though radiation has yet to decouple from matter. The evolution 
of yc is governed by the relation [50| 
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Proceeding with the mode sum as above and adopting the same approximations as above, we find that 
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where ^ls ^ 1.19 x lO"'^"^ s is the time of last scattering and Ii{m, n, a, f3, Anii„, A^ax) is once again given by Eq. (|4.35l) . 
Note that since matter-radiation equality occurs prior to last scattering, at around ^mre '^ 10^^ s, the epoch during 
which a\ decays can affect yc straddles both the RD and MD eras. Numerically evaluating the expression in Eq. (|4.40p 
from tEc to tLs, we obtain our final results for yc distortions due to late ax decay. 

In order to assess the CMB constraints on the parameter space of our bulk-axion DDM model, we now compare 
the results obtained by numerically integrating Eqs. (|4.28l) and (|4.40p with observational limits on /x and yc- The 
current limits on these quantities are [45j 



Im 



< 9 X lO"'^ , yc < 1.2 X 10" 



(4.41) 



The bound on yc for a photonic axion with c^ — 1 yields the constraints on fx, Mc and Ac shown in Fig. [6l In this 
figure, we display contours of the values of yc in {fx, Mc) space which arise in a bulk-axion model with Ac — 1 GeV 
(left panel) and with Ac = 1 TeV (right panel). In each case, we have taken £_ = gc = = 1, Hj = 1 GeV, and 
Trh — 5 MeV. Contours indicating y = Mc/mx = 1 (sohd red line) and y — {0.01,0.1,10,100} (dashed red lines) 
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FIG. 6: Contours of the CMB Conipton-j/-parameter distortion yc (black lines) produced as a result of axion decays in a 
bulk-axion DDM model with Aq = 1 GeV (left panel) and Ac = 1 TeV (right panel). In each case, we have assumed a 
photonic axion with c^ = 1 and have taken ^ — gc ~ = 1, with Hi — 1 GeV and Trh ~ 5 MeV. Contours corresponding 
to 3/ = Mc/mx ~ 1 (solid red line) and to y = {0.01, 0.1, 10, 100} (dashed red lines) are also shown. For each panel, it is 
evident that our bulk-axion DDM model amply satisfies the CMB constraints in Eq. (|4.41[) for all relevant values of fx and 
Mc, regardless of the value of y. 

have also been superimposed. For each panel, it is evident that our bulk-axion DDM model amply satisfies the CMB 
constraints for all relevant values of fx and Mc, regardless of the value of y. 

It turns out that the constraints from the corresponding bound on /i in Eq. (|4.4ip are even less stringent than those 
from the bound on yc- Thus, we conclude that both the yc-type and /i-type distortions which result from ax decays 
in our bulk-axion DDM model are well below present experimental sensitivities. Indeed, no meaningful constraint 
arises for our bulk-axion DDM model from present limits on distortions in the CMB. 

As we have discussed, neither /i nor yc can be affected by any photons which are produced by ax decays at times 
t > ^LSj after radiation and matter decouple. Such photons do, however, contribute to the diffuse photon background. 
In the next section, we will discuss the physical effects of this diffuse photon background in detail. 

G. Axion Decays and Contributions to the Diffuse X-Ray and Gamma-Ray Backgrounds 



As mentioned above, the potentially observable effects of late photoproduction from axion decays include not only 
distortions of the CMB, but also imprints on the diffuse X-ray and gamma-ray backgrounds. Observational limits 
on such imprints from instruments such as HEAO ^51,], COMPTEL [52], XMM, and Chandra [53| therefore impose 
additional constraints on the parameter space of our bulk-axion DDM model. As discussed in Sect. IIII[ there are 
two cosmological populations of decaying ax whose decays to photons can potentially leave observable imprints on 
the diffuse X-ray and gamma-ray backgrounds. The first is the population of cold axions produced by vacuum 
misalignment, which collectively compose the DDM ensemble. The second is the far smaller population of axions 
produced by interactions among the SM fields in the thermal bath after inflation. While the former population 
provides a far greater contribution to i^tot, the latter population contains a far larger proportion of heavier, more 
unstable ax, as indicated in Fig. |3l It is not clear a priori which population yields the more stringent constraint. 
Thus, it is necessary to examine the contribution to the diffuse photon background from each of these populations in 
turn. 

A photon produced at time t with initial energy E-y{t) will only contribute to the diffuse photon background if 
the universe remains transparent to electromagnetic radiation over the entire range of energies through which that 
photon redshifts as the universe evolves from t to inow A detailed analysis of the time scales and photon-energy 
ranges for which this transparency condition is attained is presented in Ref. [54J . Roughly speaking, the transparency 



window spans an energy range 1 keV < E^ < 10 TeV and a time range 10 



12 



10 



14 



< t < t„ 



with the lower limit 
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depending on the particular value of E^. Motivated by these results, we approximate the universe to be transparent 
to all photons with energies which fall within this range at all times t > iLg and opaque to all photons otherwise. 
This approximation yields a conservative bound. Moreover, we emphasize that since the dominant contribution to 
the diffuse X-ray and gamma-ray flux in our model is due to modes which decay at much later times t ^ tLS, our 
results are essentially insensitive to the precise contours chosen for the transparency window. 

The calculation of the photon flux due to late ax decays proceeds in a manner similar to the calculation of the flux 
from KK-graviton decays outlined in Ref. [2a| . The Boltzmann equation for the number density n^ of photons in the 
presence of a tower of decaying a\ takes the form 



+ 3Hn, = 2 5:BRi^-)r,^, 



(4.42) 



where once again p\ is given by Eq. ([3 



Solving this equation for n^ as a function of time, we obtain 

2£W y^ gpJ27) PA(tLs) L _ g-rA(t-tLs) 
SLS 4^ ^ A L 



(4.43) 



where s(t) is the entropy density of the universe at time i, and sls is the entropy density of the universe at the time 
of last scattering. The present-day differential energy spectrum dn-y/dEj of these photons may readily be computed 
from the relation 
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dt dz dE^ 



(4.44) 



where z is the cosmological redshift and E^ is the photon energy at redshift z. The first of these factors may be 
obtained by explicitly differentiating Eq. (j4.43p with fixed s — Snow, which yields a series of terms of the form 
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SLS 



A 



-tLs) 



(4.45) 



one for each different value of A. The second factor in Eq. (|4.44|) may be obtained by noting that the relationship 
between time and redshift during the present, matter-dominated era is well-approximated by i = inow(l + z)^^/"^ . 
Consequently, for each value of A we have 



dt 
dz 



^ ^now 



2E^ 



5/2 



(4.46) 



during the epoch of interest. The third factor in Eq. (I4.44p may be obtained by noting that each of the photons 
produced by an axion tower state ax which decays at redshift z will be monochromatic, with energy A/2, at the 
moment of decay. This implies that the present-day energies of such photons are given by -£^7(1 + z) = A/2, and hence 
that for each value of A, we have 



dz 



dE^ 



A 

2^ 



(4.47) 



Combining these expressions and summing over A, we arrive at a general formula for the contribution to the diffuse 
photon flux produced by the tower of decaying ax'- 
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Calculating the contribution to the diffuse X-ray and gamma-ray backgrounds in our bulk-axion DDM model is then 
simply a matter of applying Eq. ()4.48|) to the contribution from the two relevant populations of decaying axions 
discussed above. 

We begin by addressing the contribution from the population of axions produced by vacuum misalignment — i.e., the 
DDM ensemble itself. Once again, we focus our attention on the case of a photonic axion, for which Tx ~ r(a — > 77) 



and BR' 



(27) 



1. In this case, we find that the contribution to the present-day diffuse photon background from the 



collective decays of the ax fields is given by 
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Just as for the contributions to /i and yc in Sect. IIVFI we approximate the sum over axion modes appearing in 
Eq. (j4.49p as an integral over A. The lower limit of integration is determined by the requirement that in order for a 
photon with redshifted energy E~^ to have been produced by the decay of the axion species a\ before present day, we 
must have A > 2E^. Likewise, photons which decay before the processes which equilibrate them with the radiation 
bath freeze out will not contribute to features in the diffuse photon background. Thus, the upper limit of integration 
is set by the condition A < 2Ej{tnow /ths)'^^'^ . Furthermore, we must also require that A not exceed the cutoff scale 
/g, or be smaller than the lightest mode in the tower. Once again, we find that the resulting integral expressions can 
be written in terms of the functions Ii{m, n, a, /3, Amin, Amax) defined in Eq. (|4.35p : 
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(4.50) 



where the Af ^^ are analogous to the \™^ appearing in Eq. ()4.38p . Determining the net contribution to the differential 
photon flux from decays of the a\ for any particular choice of model parameters is thus simply a matter of numerically 
evaluating Eq. (j4.50p . 

We now turn to consider the observational limits on dn-y/dE-y. The diffuse extragalactic X-ray and gamma-ray 
background spectra have been probed by a number of experiments. In the keV — MeV region, the most current data 
are those from HEAO, COMPTEL, XMM, and Chandra; at energies above this, the most current data are those 
from EGRET and FERMI. Over this entire energy range, the diffuse photon spectrum is well-modeled by a set of 
power-law fits, and the non-observation of any discernible, sharp features in this spectrum imposes constraints on 
late relic-particle decays to photons. For the data from the COMPTEL instrument, the best power-law fit is found 
to be [521 
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while the best fit to the HEAO data is found to be 1511 
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(4.52) 



The Chandra satellite has improved upon these diffuse X-ray background constraints in the 1 keV ^ E^ < 8 keV 
range by resolving a large fraction (~ 80%)of this background into point sources. The residual spectrum in this 
region is well represented by the power law [57[ 
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In the gamma-ray region, the most stringent current limits are those from EGRET and FERMI. Data on the diffuse 
extragalactic gamma-ray background from the former instrument [55| are reliable for photon energies within the range 
1.41 GeV < .E.^ < 30 MeV, for which we find the best fit 
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(4.54) 



Note that data exist for higher photon energies as well, but given that EGRET's energy resolution is not as good at 
such high energies, and given that these data have been superseded by data from FERMI, we do not use them in 



computing this power-law fit. As for tlie FERMI data, they are well modeled by the power law [5f 
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(4.55) 



In Fig. [71 we show a set of curves (solid colored lines) depicting the total contribution to the diffuse gamma-ray 
background from the decaying a\ fields, as given in Eq. (|4.50l) . for several different values of fx within the range 
10^ — 10^^ GeV. Results are shown for Aq = 1 GeV (upper left panel), Ac = 1 TeV (upper right panel), and 
Ag = 100 TeV (lower panel). In each case, we have taken M^ = 10~" GeV, Trh = 5 MeV, and £^ = go = 9 = 1. 
In addition, we have chosen a value for Hj sufficiently large that none of the curves shown is significantly affected 
by the "inflating away" of heavy modes which begin oscillating before inflation ends. In addition to these curves. 
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FIG. 7: The diffuse photon-flux spectrum dn-,/dE-, produced from axion decays in our bulk-axion DDM model with Ag = 1 GeV 
(upper left panel), Ag ~ 1 TeV (upper right panel), and Ag ~ 100 TeV (lower panel). Each solid colored curve corresponds 
to a different choice of fx, within the range 10® — 10^® GeV. In all panels, we have taken Mc — 10^^^ GeV, Hi — 1 GeV, 
Trh = 5 MeV, and (, = go ~ 9 = i- By contrast, the dashed black contours represent the upper bounds on dn^/dE^ derived 
from observational limits on the diffuse photon flux using a number of instruments sensitive in the X-ray and gamma-ray 
regions. As evident from these plots, the diffuse-photon-background contribution arising from axion decay in our bulk-axion 
DDM model is consistent with all observational limits when Ag is large. 
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we also display contours corresponding to the upper limits on the diffuse X-ray and gamma-ray fluxes (black dashed 
lines) given in Eqs. (|4.5ip through (|4.55l) . Any choice of model parameters for which the differential photon flux 
dn-y/dE^ exceeds any one of these observational-limit contours for any value of E^ is excluded. The results shown 
in Fig. [7] indicate that while it is not trivial to satisfy these observational limits in our bulk-axion DDM model, the 
contributions to the diffuse X-ray and gamma-ray fluxes from a\ decay are indeed sufficiently small that these limits 
are satisfied when Ag is large. 

We now consider the contribution to drij/dEj from the population of axions generated by their interactions with 
SM fields in the thermal bath after inflation. The contribution to the diffuse photon flux spectrum drij/dEj generated 
by such a population of axions is once again given by Eq. (j4.48p . but with p\(ti^s) in Eq. (13. 8p now replaced by 
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(4.56) 



as follows from Eq. (J3.42I) . To derive an estimate for the expected contribution to driy/dE-y from the resulting equation, 
we proceed in essentially the same way as we did in calculating the contribution from axions produced via vacuum 
misalignment. The results of this calculation are shown in Fig. |H]for parameter values within or near the preferred 
region of parameter space for our bulk-axion DDM model. Speciflcally, we have taken Mc = 10^^^ GeV, Aq = 1 TeV, 
Irh = 5 MeV, and ^ = gc = 1. The solid colored curves shown correspond to several different choices of fx ranging 
from fx = 10^^ GeV to fx = 10^^ GeV. Once again, the dashed black lines indicate the observational limits on 
additional contributions to dn^/dE^. It is clear from Fig.[5]that while the contribution to the diffuse X-ray flux from 
thermal axions within our preferred region of parameter space is certainly not negligible, it is also consistent with 
current observational limits. We therefore conclude that even after the contribution from thermal axions is included, 
our bulk-axion DDM model is consistent with X-ray and gamma-ray data. 
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FIG. 8: The diffuse photon flux spectrum dn^/dE-, produced from the decays of a population of ax produced by interactions 
among the SM particles in the thermal bath after inflation. The left panel shows the results for Ag ~ 1 TeV, while the right 
panel shows the results for Ag = 100 TeV. In each case, we have taken Ale = 10~^^ GeV, Trh ~ 5 MeV, and ^ = gtG = 1. 
The solid colored curves indicate the diffuse-photon-flux contributions corresponding to different choices of fx ■ As in Fig. [71 
the dashed black contours indicate the upper bounds on dn-y/dE^ derived from observational limits on the diffuse X-ray and 
gamma-ray fluxes, and we see that our model is consistent with these bounds. 



H. Axion Decays and Big-Bang Nucleosynthesis 



The accord between the primordial abundances of light nuclei inferred from observation and the predictions for those 
abundances within the framework of standard BBN has been one of the greatest triumphs of theoretical cosmology. 
However, these predictions depend sensitively on the cosmological parameters during the nucleosynthesis epoch. For 
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example, the presence of additional relativistic degrees of freedom in the thermal bath during BBN can substantially 
distort the abundances of the light elements away from their observed values. In addition, the decays of unstable 
particles during or after the BBN epoch can also alter these abundances via the injection of both entropy and energy 
into the thermal bath. We must therefore ensure that the collective effects of a\ decays in our model are sufficiently 
small so as not to disrupt the successful generation of light-element abundances via standard BBN. 

Limits on the abundance of a single unstable relic particle x from BBN are typically phrased as bounds on the 
number density n* that x would have at present time if it were absolutely stable. In general, the BBN bound on 
n* for any given relic particle depends on the lifetime t^ of that particle. The most stringent limits are obtained for 
lifetimes r^ ~ 0(10^ — 10^° s), for which the corresponding constraint is roughly |58ll59l| 
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where 
form 



410.5 cm ^ denotes the present-day number density of photons. This limit can also be written in the 



fi* < 1.7 X 10"^ 



(4.58) 



where il* is the relic abundance that x would have at present time if it were absolutely stable. 

Once again, however, as with other constraints on traditional models of decaying dark matter (such as those from 
the CMB and the diffuse X-ray and gamma-ray backgrounds), these constraints are not readily applicable to models 
within the DDM paradigm, since the dark-matter candidate in these models is an ensemble with no single, well-defined 
mass or lifetime. Thus, we must reexamine the derivation of the BBN constraints on decaying relic particles in order 
to determine what restrictions these considerations place on the parameter space of our bulk-axion DDM model. 
While a detailed calculation of the precise limits BBN considerations impose on DDM scenarios in general is beyond 
the scope of this paper, it is straightforward to demonstrate that BBN constraints do not significantly restrict the 
parameter space of the particular model which concerns us here. 

We begin by noting that in traditional, single-particle dark-matter scenarios, an unstable dark-matter candidate 
X with a relic abundance il^ ~ ^cdm is generally consistent with all astrophysical and cosmological limits on dark- 
matter decays, provided that Ty. > r™™ ~ 10^^ s [53|. It therefore follows that any ax in the DDM ensemble with 
a lifetime t\ > r™'" will have no impact on BBN within regions of parameter space in which the WMAP constraint 
f^tot ^ f^CDM on the total dark-matter relic abundance is satisfied. Thus, we may safely conclude that our bulk-axion 
model of dynamical dark matter is consistent with BBN constraints within such regions of parameter space, provided 
that 
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FIG. 9: Contours of the collective contribution to Q*^^ from the set of ax with lifetimes tx < t™'" for a DDM ensemble 
of photonic axions with Cy = 1. The left, center, and right panels display the results for Ag ~ 1 GeV, Ag = 1 TeV, and 
Ag = 100 TeV, respectively. In each case, we have taken Trh ~ 5 MeV, Hi = 100 TeV, and ^ = go = d = 1. In each panel, 
we see that BBN constraints are amply satisfied throughout essentially the entire region of parameter space shown. 
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where ^*^^ denotes the collective contribution which the set of a\ with lifetimes t\ < r™'" would have made to the 
dark-matter relic abundance at present time if they were absolutely stable. In other words, the BBN constraint we are 
imposing in Eq. (|4.59|) effectively rests upon the extremely conservative approach of treating all states in the DDM 
ensemble whose lifetimes are less than t™'" as if they had lifetimes tx which are in the range which is most dangerous 
for BBN, namely t\ ^ 10^ — 10^° s. We emphasize that while this criterion is a sufficient condition for successful 
BBN, it does not represent the true BBN constraint, which is always far less stringent. 

In Fig. [9l we display contours of Cll^^. for a DDM ensemble of photonic axions with c-y = 1, as a function of fx, Mc, 
and Aq. The left panel shows the results for Aq — 1 GeV, the center panel for Aq — 1 TeV, and the right panel for 
Ag — 100 TcV. In each case, we have taken Trh = 5 MeV, Hj = 100 TeV, and £^ = go = 9 = 1. In each panel of 
Fig. [HI we see that the criterion in Eq. (j4.59|) is amply satisfied throughout essentially the entire region of parameter 
space shown. It therefore follows that our bulk axion model is consistent with successful BBN throughout this region 
of parameter space. 

I. Axion Decays and Late Entropy Production 

One additional physical consequence of the late decays of unstable relic particles is the generation of entropy as those 
particles "dump" their energy density into the radiation bath. Indeed, a number of considerations place constraints 
on late entropy production from decaying particles. For example, late entropy generation can upset the light-element 
predictions from standard BBN and produce observable features in the CMB. In this section, we examine the effect 
of the late decays of the ax on the entropy density of the universe in our bulk-axion DDM model as a function of time 
in order to verify that no perceptible effects can arise which might serve to exclude our model. 

During any given epoch, the entropy density of the universe is dominated by the contribution from radiation and 
therefore well approximated by 

^- E^ = ^9.s{T)T\ (4.60) 

i 

where the index i runs over all relativistic particle species, Ti is the temperature associated with any particular such 
species, and g*s is the number of interacting degrees of freedom at temperature T. During the early stages of the 
history of the universe (prior to neutrino decoupling), all such species are characterized by a common temperature 
Ti « T. During such epochs, g*s{T) « g*{T), and the entropy density is therefore directly proportional to the total 
energy density Prad of radiation. Indeed, even during subsequent epochs, (?*s and g* remain roughly similar, and prad 
remains a good indicator of the entropy density. Thus, by evaluating the contribution to prad from ax decays in our 
bulk axion DDM model, we can assess the effect of these decays on both the energy and entropy densities of the 
universe. 

In the LTR cosmology, as in the standard cosmology, Prad evolves according to an equation similar to Eq. (|4.30p : 

^^rad _ .„ , r ^ _L Y^ RTJ (fad) 



dt 



iHp,,d + r^p^ + J2 ^^r^'^xPx ■ (4.61) 



This equation assumes the presence of a tower of decaying oa, where BRj^ is the total branching fraction of ax 
into relativistic particles. Note, however, that since we are working within the context of LTR cosmology, the effects 
of inflaton decays on the energy and entropy densities of the universe remain relevant until very late times t ^ tRH- 
Thus we have explicitly included an additional source term T^p^ in Eq. (J4.61I) to account for the effect of such inflaton 
decays, where F^ and p^ respectively denote the decay rate and energy density of the inflaton field (f). 

The contribution to prad from inflaton decays can readily be calculated from standard results pertaining to the LTR 
cosmology (for a review, see, e.g., Ref. [lal)- As the universe exits the inflationary epoch at a time tj w 2/(3iJ/), the 
energy density stored in the inflaton field is initially p^ = Pdit — 3HJMp. During subsequent epochs, the inflaton 
source term for radiation is approximately given by 

where k, is defined as in Eq. p.3p . and we have used the fact that the inflaton-decay rate is related to the reheating 
time by F^ « l/{2t^^). Note that this source term is negligible at times t ^ %,Hi when by definition p^ <^ Pradi and 
hence can safely be neglected at such times. By contrast, at early times t < ^rh, the inflaton source term is expected 
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to dominate in Eq. (|4.61|) . in the sense that 



r^P0 » J^BRf'^^rAPA 



(4.63) 



Whenever this condition is satisfied, the contribution to prad from a\ decays is inconsequential compared to that from 
inflaton decays, and the axion source term can therefore safely be neglected. 

In assessing the contribution from ax decay, we once again choose to focus on the case of a photonic axion with 
c-y — I; this implies that the decay mode ax — > 77 dominates the contribution to prad- In this case, the source term 
for radiation due to a\ decay is just the source term for photons given in Eq. (|4.37p . In this case, solving Eq. (I4.6ip 
for prad, we find that 
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(4.64) 



where Prad(0 i^ the solution for Prad(^) in the absence of any additional contribution from ax decays. Once again 
making use of the integral functions /^(to,^, a,/3, Amin, Amax) defined in Eq. (|4.35p to approximate the sum over 
modes, we obtain 
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In Fig. I10[ we show how the contribution to p^ad from ax decays in our bulk-axion DDM model evolves with time 
for a variety of different choices of model parameters. The left panel shows results for Aq — 1 GeV, the center panel 
for Ag = 1 TeV, and the right panel for Aq — 100 TeV. The solid colored curves in each panel correspond to different 
choices of fx within the range 10^° — 10^® GeV. For all curves shown, we have assumed a photonic axion with c^ = 1, 
and we have taken M^ = 10"" GeV, Trh = 5 MeV, Hj = 100 TeV, and ^ = go = = 1. The black dashed curve 
represents the total value of Pradj which includes the standard contribution from inflaton decays during the reheating 
epoch. Since such inflaton decays constitute the dominant source for radiation prior to the end of reheating, the 
range of times shown in each panel extends from i^H to present time. The value of Hj has been chosen here to be 
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FIG. 10: The total contribution to the radiation energy density prad from photonic a\ decays in our bulk-axion DDM model 
(solid lines) , plotted as functions of time for a variety of different choices of fx ■ The left panel shows the results for Ag = 1 GeV, 
the center panel for Ag = 1 TeV, and the right panel for Ag = 100 TeV. In each case, we have assumed a photonic axion 
with ^ = gc = 9 ~ 1, and we have taken Mc = 10~^^ GeV, Trh = 5 MeV, and Hi = 100 TeV. Also shown in each panel is 
the total value of prad as a function of time in the LTR cosmology (black dashed line), which includes the contribution from 
inflaton decay. In all cases, the collective contribution to prad from ax decays at all times t < inow remains negligible compared 
to the primordial contribution generated via inflaton decays during reheating. Thus our bulk-axion DDM model does not lead 
to overproduction of either radiation-energy density or entropy during any prior cosmological epoch. 
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sufficiently large that the effect of heavier a\ with A > 3Hj/2 being inflated away is unimportant. Note, however, 
that for significantly smaller values of Hj, the contribution to p,ad from axion decays can be further suppressed by 
this effect. 

The differences among the curves shown in Fig. [TU] for different choices of fx and Aq ultimately stem from the 
effects of axion mixing on the abundances px and decay widths Tx of the individual axion modes. The results shown 
in the left panel correspond to the case in which Aq is sufficiently small that y ^ 1 for all choices of fx shown. 
In this small-mixing regime, A > Atrans for all but the lowest-lying mode in the axion KK tower, and Eqs. ()2.12p 
and p.Sp imply that px oc f^^ and Tx oc f^^. It therefore follows that the photon source term BRj^ ^ Txpx associated 
with each ax within this regime decreases uniformly and substantially with increasing fx, as indicated. By contrast, 
as Ag is increased, several competing effects play an increasingly important role in determining the magnitude of 
BR)^ '^ Txpx for certain A. This is because Atrans increases with increasing Ag; hence for large Ag a greater number of 
the ax are brought into the A < Atrans regime, in which px oc fx and Tx (x fx- Increasing fx therefore has the effect 
of increasing the initial magnitude of the photon source terms associated with the ax in this regime. However, the 
lifetimes of these modes also increase with increasing fx , and hence the transfer of their energy density to radiation is 
deferred until later times, when p^ad is smaller and the contribution from ax decays can have a proportionally greater 
impact. The interplay between these effects results in the behavior shown in the right two panels of Fig. [TU] 

Note that the curves for the total energy density shown in Fig. 1101 which are dominated by the contribution from 
inflaton dynamics, drop more rapidly as a function of time than the contributions from axion dynamics. This reflects 
the continuing generation of new radiation energy density from the ongoing decays of the individual ax within our 
DDM ensemble. In all cases, however, the collective contribution to prad from ax decays at all times t < t^ow remains 
negligible compared to the primordial contribution generated via inflaton decays during reheating. Thus our bulk- 
axion DDM model does not lead to overproduction of either radiation-energy density or entropy during any prior 
cosmological epoch. 

J. Vacuum Energy and Overclosure 

In traditional dark-matter scenarios involving a single, stable dark-matter candidate x, the dark-matter relic abun- 
dance n,^ increases monotonically up to and beyond the present time. As a result, verifying that ily^ satisfies WMAP 
constraints at the present time is sufficient to guarantee that x does not overdose (or prematurely matter-dominate) 
the universe at all previous times as well. However, one of the hallmarks of the DDM scenario is that this is no 
longer true: although iltot likewise experiences a Hubble-driven growth during the earliest phases of the evolution of 
the universe, this quantity can nevertheless drop during later epochs. This is possible within the DDM framework 
because the single, stable dark-matter candidate x characteristic of most traditional dark-matter scenarios is replaced 
by a complex, multi-component dark-matter ensemble whose constituents can have a broad spectrum of lifetimes and 
abundances. As a result, the decays of certain dark-matter components within the ensemble can cause fltot to decline 
— even prior to the present day. Indeed, such behavior for iltot can be quite dramatic, and is illustrated in Fig. 6 
of Ref. [l[ for the special case in which the DDM ensemble consists of a KK tower of decaying dark fields. Thus, 
within the DDM framework, it is no longer sufficient to verify that iltot satisfies overclosure constraints at the present 
time; we must also verify that it has satisfied such overclosure constraints (and constraints from premature matter- 
or vacuum-energy domination) at all prior moments during the history of the universe. 

It turns out, however, that this is not a problem in our bulk-axion DDM model. Since our model already satisfies 
WMAP constraints at present time within our preferred region of parameter space [4I , it can run afoul of overclosure 
constraints in the past only if the negative rate of change of fitot is sufficientlygreat that fltot rnight have exceeded 
unity within the past history of the universe. However, as discussed in Refs. [l|, yl; this rate of change is described 
by an effective equation-of-state parameter Wcs, and two things are already known about the value of this parameter 
in our model: first, it is extremely small at the present day, i.e., 10~^^ < WcS ^ 10~^^ "^l, and second, it was even 
smaller in the past. Indeed, this latter assertion follows from the generic behavior of WcS shown in Fig. 8 of Ref. [1|: 
for a generic KK tower, WcS reaches its maximum at the present day and is exponentially smaller prior to this time. 
Thus, working backwards from the present epoch, and given the finite age of the universe, we see that it is not possible 
for J7tot to have violated overclosure bounds at any point during the history of the universe. 

One related concern which arises in our bulk-axion DDM model, due to our reliance on the misalignment mechanism 
for the generation of the primordial relic abundances of the ax is the risk of premature vacuum domination. Indeed, 
any ax for which tx > tQ will contribute to the total dark-energy abundance fivac during the period when Iq ^t<tx, 
within which its energy density px is non-vanishing but before which it begins oscillating. Since px remains constant 
during this period, the contribution to Ovac scales like fix (x t'^ during any MD or RD epoch. Since this represents a 
rate of increase far faster than that associated with matter or radiation, the threat of premature vacuum domination 
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from fields which remain as vacuum energy for a long duration is of particular concern. Indeed, in extreme cases, 
such fields could potentially give rise to an additional period of inflation, leading to gross inconsistencies with the 
predictions of BBN, CMB data, and so forth. 

In our bulk-axion model, however, it is straightforward to demonstrate that no such inconsistencies with obser- 
vational data arise. The masses of all of the ax, with the sole exception of the zero mode oq, are bounded from 
below by the Newton's-law-modification constraint in Eq. (|4.ip . since Xi > AIc/2 for i > 0. For all such modes with 
^A > to, this constraint on A implies a bound t\ > 6.75 x 10"-'^'* s on the oscillation-onset time of the mode. (The 
remaining modes, for which tc = t\, never contribute to r^vac-) This time scale is sufficiently early that the collective 
vacuum-energy contribution from these ax poses no threat of overclosure or premature vacuum-domination. The fix 
contributions from these fields simply do not have time to grow to a problematic size. 

This leaves only the contribution from ag, whose oscillation time scale can be substantially longer than the upper 
limit quoted above for the higher modes in situations in which y ^ 1. Since Ax^ ~ 1 in this limit, Eq. p.8p implies 
that prior to the time txg at which it begins oscillating, the relic abundance of oo is given by 

^^« " 2^ (;:) ■ ^'-''^ 

Therefore, one finds that by the time of oscillation, which is given by tx^ ~ Kxo/2mx in this limit, flx^ will have 
grown to 
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8 A/2 



^Ao(iAo) « 0T7T • (4-67) 



This result is independent of mx , and implies that the contribution of the ao to fivac is not a cause for concern for sub- 
Planckian values of fx- Indeed, this is to be expected: in this regime, ao functions effectively like a four-dimensional 
axion. Early vacuum-energy domination is known not to be a problem for light axions and axion-like particles (see 
Ref. [35| and references therein) in purely four-dimensional theories. 

K. Misalignment Production and Isocurvature Perturbations 

In an inflationary cosmology, fluctuations in the energy density of any population of particles produced thermally, 
i.e., via rapid interactions in the radiation bath during the reheating phase, stem from the primordial perturbations in 
the energy density of the inflaton field. Consequently, such fluctuations are of the so-called adiabatic type — that is, 
they represent spatial variations in the total energy density, but not in the relative contributions of individual particle 
species to that total density. Such variations, in turn, imply fluctuations in the local spacetime curvature and are 
therefore sometimes also referred to as curvature perturbations. By contrast, fluctuations in the energy density of any 
population of particles produced via means uncorrelated with the inflaton field (and therefore non-thermal) can also 
give rise to fluctuations of the isocurvature type — i.e., perturbations in the relative contributions of different species 
to the total energy density, with that total energy density held fixed. Recent WMAP observations of the CMB power 
spectrum, taken in combination with baryon acoustic oscillation (BAO) measurements and supernova data, place a 
stringent bound [lO[ on any deviations from adiabaticity in primordial energy-density fluctuations. This bound is 
typically expressed in terms of the fractional contribution ao to the CMB power spectrum from axion isocurvature 
perturbations: 

"° ^ nZZ^}i < 0-0^2 , (4.68) 

where {{5T/T)1^^) and {{5T /T)f^^ respectively denote the total average root-mean-squared fluctuation in the CMB 
temperature, and the average root-mean-squared temperature fluctuation due to isocurvature perturbations alone. 
Since the ax fields which compose our dynamical dark-matter ensemble are presumed to be produced non-thermally, 
via the misalignment mechanism, it is necessary to investigate the implications of this bound for our model. 

Our discussion of isocurvature perturbations in our bulk-axion DDM model in large part parallels the discussion 
of such perturbations in traditional QCD axion models presented in Ref. [60], to which we refer the reader for a 
more complete introduction and discussion of the formalism and methodologies used. It turns out to be convenient 
to express the fluctuations of any given ax in terms of the fractional change Sx in the ratio of its number density nx 
to the entropy density s of the universe. This quantity can be written in the form 

5{nxls) dnx ST 

Sx ^ -r^ = ^ - 3— . 4.69 

(nx/s) nx T 
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We assume that the production of all other particle species i/ji (i-C, the SM fields) ultimately results from inflaton 
decay, and that the density fluctuations for these species are purely adiabatic, with Si — 0. Since, by definition, the 
fluctuation 5p in the total energy density vanishes for isocurvature fluctuations, it therefore follows that the sum of 
the fluctuations in the energy densities of the various particle species obeys a constraint which may be written in the 
form 

J2p,(^S,+3^-^^+3Y,p~+4p,J-f = 0, (4.70) 

where the pi denote the energy densities associated with massive species other than the oa, and prad once again 
denotes the total energy density of radiation. In our bulk-axion DDM model, the abundances of all of the a\ are 
determined by a single misalignment angle 6. As discussed in Ref. [l|, this reflects the ultimate five-dimensional 
nature of the axion field. This in turn implies that the density fluctuations Snx for all of these fields are determined 
by the fiuctuations 66 in this misalignment angle generated by quantum fiuctuations during inflation. The fact that 
the fluctuations 5n\ are all determined by 69 implies that the S\ = S are essentially equal for all ay, hence Eq. (|4.70|) 
simplifies to 

^totS = - 3 ('anat + ^n,,d) ^ , (4.71) 

where fimat denotes the total abundance of matter in the universe, including the contributions from baryonic matter, 
the ensemble of dark axions, and any other particles which might contribute to the dark-matter relic abundance, 
and r^rad is the relic-abundance contribution from radiation. This expression is identical to that which describes the 
isocurvature perturbations associated with a single, four-dimensional axion. Therefore, assuming that the fiuctuations 
in 9 are Gaussian, it follows that in our axion DDM model, ag is given by the standard expression [60| 

8 fn;^^Y 1 a^{29^ + a^) 

"°-25U;;.aty mmi,) {9^+aiY ' ^^ ^ 

where O^^t denotes the present-day value of fimat, and where Cg = {{56)'^) denotes the variance associated with 
fluctuations in 9. 

This result makes intuitive sense. Although our DDM model has essentially partitioned the total dark-matter 
abundance amongst a large number of different KK axion fields, the underlying flve-dimensional nature of the KK 
tower has correlated the individual fluctuations of these flelds so that they are governed by the fluctuation of a single 
misalignment angle 9. It is therefore not a surprise that the expected magnitude for isocurvature fluctuations in our 
model turns out to be no greater than it is standard, four-dimensional axion models. 

All that remains, then, for us to do in order to determine the value of ag in our bulk-axion DDM model, is to assess 
the magnitude of a^ . Assuming again that the fluctuations in 9 are Gaussian, this quantity is given by 

-I = -It . (4.73) 

Since we are operating within the context of an LTR cosmology with Trh ~ 0{MeV), as discussed above, it is by no 
means problematic (and in fact quite natural) for Hj ^ fx- Therefore, as long as 6* '-^ C'(l), as might be expected 
from naturalness considerations, it can safely be assumed that 9 ^ ae- Substituting into Eq. (I4.72[) the experimentally 
observed 10] values {{5T/T)1^^) sa (1.1 x 10^^)^ and r^J^^^t ~ 0.262 we flnd that ao is weU approximated by 

ao « 1.95 X 10^ (^^■^B^\ (4.74) 

in our bulk-axion model. Combining this result with the upper bound on ao quoted in Eq. (I4.68|) yields the constraint 

Hi < 6.07 X 10~^ ( Ur^] ■ (4.75) 

We consider the case in which fl^^^ ~ ^cdm and in which the axion ensemble is responsible for essentially the entirety 
of the observed dark-matter relic abundance. This ocrresponds to fx ~ 10^"* — 10^^ GeV. We then find that for 
9 ^ 0(1), the resulting constraint Hj < 10^ — 10^° GeV on the Hubble parameter during infiation is relatively mild. 
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Indeed, there is no difficulty in satisfying tfiis constraint in eitfier tlie standard or tfie LTR cosmology. We thus 
conclude that isocurvature perturbations do not present any problem for our bulk-axion model of dynamical dark 
matter. Moreover, a low scale for Hj can be regarded as natural in the context of an LTR cosmology. 

It is worth remarking, however, that the above results have implications for the detection of primordial gravitational 
waves. Limits on primordial gravitational waves from observations of the CMB can be conveniently parametrized in 
terms of the scalar-to-tensor ratio r. For example, consider single-field models of inflation, in which r — 16e, where 
e = Mp{V' /V)^ / {4:Tt) is the inflaton slow-roll parameter, with V and V denoting the inflaton potential and its first 
derivative with respect to the inflaton field, respectively [61]. In the context of our bulk-axion DDM model, the 
standard relation (see, e.g., Ref. jlO| ) between r and ag takes the form 



Ml I n* J 1 - ao 



(4.76) 



As discussed above, consistency with the bounds in Eqs. (|4.68l) and (I4.75|) requires that Hj <C 2Tr fxO and ao <C 1. In 
this regime, one finds that the expected tensor-to-scalar ratio is essentially independent of ^l^t ^^d well approximated 

by 

Current WMAP observations, again in conjunction from BAO and supernova data, place an upper bound r < 0.22 
on the tensor-to-scalar ratio [1^. Thus, Eq. (|4.77p results in a constraint Hj < 6.7 x lO^'^ GeV on the Hubble scale 
during inflation — a constraint which Eq. ()4.75|) implies is already automatically satisfied, even for 0(1) values of the 
misalignment angle 6. The upshot is therefore that while there is no conflict between current limits on isocurvature 
perturbations and the predictions of our bulk-axion DDM model, the requirement that Hj be relatively small in this 
model suggests that r should likewise be quite small — at least in the simplest of inflationary scenarios. Constraints 
on the spectral index Ug from WMAP [IQ\ can simultaneously be satisfied for small r without difficulty, for example 
in negative-curvature models of inflation, which tend to predict small r |62| . 

In summary, we conclude that current constraints on isocurvature perturbations can be satisfied in our bulk-axion 
DDM model without too much difficulty. However, we note that any conclusive measurement of r within the sensitivity 
range of the Planck satellite would have severe ramifications for this model. 

L. Axion Abundances and Quantum Fluctuations During Inflation 

Thus far in this paper, we have disregarded the effects of the quantum fluctuations that naturally arise for any 
massless or nearly massless field during the infiationary epoch. In particular, the low-momentum modes of any a\ 
in our model with a mass A < Hj have wavelengths which exceed the Hubble length during inflation; excitations of 
such low-momentum modes are therefore indistinguishable from a VEV and consequently do not inflate away. These 
excitations necessarily yield a primordial energy-density contribution in our bulk-axion DDM model which cannot 
be avoided in any inflationary cosmology. Consistency with the relic-abundance predictions discussed in Sect. IIII Al 
therefore requires that this primordial energy density be small compared to that which results from misalignment 
production. 

In particular, it is possible to formulate a condition that ensures that these quantum fluctuations not invalidate 
our previous analysis. Clearly, one criterion that any such condition must enforce is that such fluctuations not have a 
significant effect on the total relic abundance of the ensemble. We may formulate this constraint as a requirement that 
the difference between the full present-day relic abundance ri^^^, which incorporates the effect of these fiuctuations, 
and the result il'^^^. obtained in the absence of such corrections be negligible — i.e., that 

|f^:ot-^t*ot| « ^tof (4.78) 

While the condition in Eq. (|4.78|) is certainly a necessary one, it is not by itself sufficient to ensure that vacuum 
fluctuations during inflation do not lead to phenomenological difficulties for our model. This is because within DDM 
framework, dark- matter stability is not a requirement, and consistency with observational constraints is arranged by 
balancing decay widths against abundances across the entire dark-matter ensemble. Indeed, as we have demonstrated, 
misalignment production provides precisely the right relationship between the Q\ and Tx to mitigate the deleterious 
effects of the heavier, more unstable states in our ensemble and render our model phenomenologically viable. We 
must therefore ensure that this delicate balance is not disrupted by the effects of vacuum fluctuations during inflation. 
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Within the preferred region of parameter space of our bulk-axion DDM model, as discussed in Sect. IIII Al the 
osciUation-onset times for the hghter a\ in the tower are staggered in time. As a result, these lighter modes collectively 
dominate in Jltot- It therefore follows that whether or not the total- relic-abundance constraint in Eq. (|4.78|) is satisfied 
depends primarily on how vacuum fluctuations affect the abundances of these most abundant modes alone. By 
contrast, the balancing of lifetimes against abundances depends on the properties of the full KK tower, and not 
merely on the attributes of the lighter modes which dominate J^tot • The corresponding condition we impose on our 
model therefore represents an even stronger constraint than the one appearing in Eq. (|4.78l) and indeed subsumes it. 
To wit, we require that the full relic abundance fix of each axion mode not differ significantly from the corresponding 
abundance Qx obtained in the absence of corrections due to vacuum fluctuations during inflation — i.e., that 

|f^A-f^A| < ^\ for aU A . (4.79) 

We emphasize that this is an overly conservative constraint, and that consistency with observational data is certainly 
possible even if vacuum fluctuations do have a signiflcant effect on the abundances of certain ax. However, as we shall 
demonstrate, the restriction that this overly conservative constraint imposes on our model (which primarily turns out 
to take the form of an upper bound on Hj for any allowed choice of fx, Mc, and Aq) is not terribly severe. 

In order to determine how this condition restricts the parameter space of our model, we must first assess what 
effect vacuum fiuctuations during inflation have on the individual energy densities px and relic abundances fix of the 
constituent fields in our dark-matter ensemble. We begin by noting a generic result in inflationary cosmologies (for 
a review, see Ref. '63']), namely that the variance ((/>^) in the amplitude of any light scalar (p with a mass m^ < Hi 
induced by vacuum fluctuations during inflation is given by 

(0^) - ^ , (4.80) 

where Sti denotes the duration of inflation. A fluctuation of this order will therefore be induced in the amplitude 
of any axion in our dark-matter ensemble with a mass smaller than Hj. Moreover, we note that the relationship 
between Stj and Hj is constrained by the fact that successful resolution of the smoothness and flatness problems 
requires N^, « HiStj > 60, where Ne denotes the number of e-foldings of inflation. In typical scenarios, N^ lies only 
slightly above this lower bound; hence Stj is typically expected to be such that HiStj ^ 0(100). Wc will frequently 
express our results in terms of N^ in what follows. 

We begin our discussion the effect of these fluctuations on the abundances of the constituent particles in our dark- 
matter ensemble by examining the simple case in which tc ^ tj. In this case, the axion mass-squared matrix attains 
its asymptotic, late-time form before inflation ends, and the ax are consequently already the axion mass eigenstates 
during the inflationary epoch. Thus, we find that the total energy density associated with each ax with A < Hj at 
the end of inflation is given by 

P^(tI) ^ lx'(0Ax.fx+vx^^4^) , (4.81) 



where 



2 



vx - <' ' , :: ,; (4.82) 




is a random coefficient of which parametrizes the fluctuation in the field ax . 

Before proceeding further, we remark that the above results depend critically on the assumption that Iq ^ tj. In 
other words, we have assumed that the instanton dynamics associated with the gauge group G has already occurred 
and made its contributions to the KK masses prior to the onset of the quantum fluctuations that arise due to inflation. 
By contrast, if to > i/, the quantum fluctuations will occur flrst, when the axion mass matrix is still diagonal and 
when the KK momentum modes and mass eigenstates coincide. In such cases, these are the modes which develop 
quantum fluctuations, and the mode-mixing induced by the instanton dynamics occurs only later. 

This distinction is important, because the resulting energy density for each ax takes a somewhat different form 
when tQ ^tj: 
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Pa = ttA 
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y^ Ux7i I OfxSnfi +rin 



27r 



(4.83) 



In this expression, Uxn is the mixing matrix in Eq. (|2.14p and ri„ is the analogue of rjx discussed above, with T]n taking 
non-zero values only when n < Hj/Mc. 
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A priori, this expression results in a different value for p\ than that in Eq. (I4.81[) . However, it turns out that the 
eventual constraints associated with Eq. (|4.83l) are no more stringent than those which we shall eventually calculate 
for Eq. ()4.81|) . In order to understand why this is the case, let us consider an even more dramatic situation in which 
r]n actually takes a fixed, positive value rj for all n — even values of n beyond the infiationary cutoff Hi/Mc- In this 
case, we can make use of the identity 



J2 Uxn = ./(A) A, 



(4.84) 



ra=0 



where /(A) = (A^ + \/2 - l)/\/2, in order to rewrite Eq. (|483l) in the form 
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:X' 



eAjx + f{X)v"'^^' 



2tt 



(4.85) 



Remarkably, this is essentially the same expression as we would have obtained from Eq. (|4.8ip when rjx — rj for all 
A, except that the fluctuation contribution now comes multiplied by an extra "scaling" factor /(A). It is easy to 
verify that /(A) — ?► 1 as A — >■ oo, whereas for small A we find that /(A) ^ 1. This indicates that the effects of the 
inflation-related quantum fluctuations are actually suppressed for the lighter modes, relative to what occurs in the 
case with to ^ tj. The magnitude of this suppression depends on y, and is more severe when y ^ I {i.e., when the 
axion modes are more fully mixed). We thus conclude that the contributions from the quantum fluctuations that arise 
during inflation are greater when they occur after the instanton dynamics turns on (and after the KK mode-mixing) , 
rather than before. We shall therefore concentrate on the ta ^ tj case in what follows. 

Given the result in Eq. (|4.8ip , we see that the effect of vacuum fluctuations on the p\ will be small for values of A 
which satisfy the condition 



OAxh 
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HiVK 

2tt 



(4.86) 



Since Ax is a monotonically decreasing function of A, it follows that within any given tower of ax, there exists a 
critical mass value 



Aflu 
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V2 




NeHJ 



y^ 



1/2 



(4.87) 



below which the effect of vacuum fluctuations on the corresponding energy density px is negligible. These px are 
therefore well approximated by Eq. (|3.8|) . and the corresponding abundances fix are given by Eq. (|3.9p or Eq. (|3.10p . 
depending on the value of ^a- By contrast, for Afluc ^ A < Hj, the effect of vacuum fluctuations overwhelms the effect 
of vacuum misalignment. The initial energy density of each ax in this regime is therefore effectively set at tj and is 
approximately given by 



Pxitj) « ^X'H] 



(4.88) 

Since the Newton's-law-modiflcation bound on Mc in Eq. (|4.ip implies that tx ^ inn for each such field, it therefore 
follows that at all subsequent times, the corresponding relic abundance is given by 

1/A < t < tRH 

O, ^ '^;'r.\i e'''^^'-''^ X <! 7T 1 ^ ) tRH < i < %RE (4.89) 

^ 4 \ tRH / 

To summarize, we see that the axion KK tower separates into three distinct regimes within each of which different 
physics plays a principal role in determining 17 a- In the A < Afluc regime, the effect of vacuum fluctuations on J7a 
is negligible and the results in Sect. IIII Al continue to hold. In the Afluc ^ A < Hj regime, the opposite is true: 
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FIG. 11: Contours of the critical value _ff"'' in {fx,Mc) parameter space. As discussed in the text, choosing Hi <C i/"'* 
guarantees that misalignment production dominates over vacuum fluctuations in determining the relic abundance fi^ of all ax 
in our DDM ensemble, as desired. Here, we have taken Tan = 5 MeV, A^e = 100, and £, — go = ~ ^ and assumed a photonic 
axion with c^ = 1. The left, center, and right panels display the results for Ag = 1 GeV, Ag = 1 TeV, and Ac = 100 TeV, 
respectively. 



vacuum fluctuations dominate and the abundances of the a\ are given by Eq. (|4.89l) . Finally, in the \> Hi regime, 
the wavelengths of even the lowest-lying momentum modes of each a\ fall short of the Hubble length during the 
inflationary epoch. Such modes therefore behave unambiguously like particles, and are consequently inflated away. 

We are now ready to address the constraint we have imposed on the individual abundances I^a in Eq. (14.791) . Since 
the effect of vacuum fluctuations is negligible both for A > Hi and for A < Afluc, it follows that this constraint 
will be satisfied whenever Hi <^ Afluc. Moreover, since Afluc itself decreases with increasing Hi, as indicated in 
Eq. (|4.87p . we find that our constraint may be expressed in the form Hi <^ i/"'', where Hf"^ is the value of the 
Hubble parameter during inflation for which Hi — Afluc. In Fig. [TlJ we display contours of Hf"^ as a function of the 
model parameters fx, Mc, and Aq. For the large values of Ac characteristic of our preferred region of parameter 
space, we observe that the constraint in Eq. (|4.79|) is satisfied for Hi ^ i/|"' ^ 0(10— 100 GeV). For smaller values 
of Ag, although the constraint is certainly more severe, we nevertheless observe that the bound can be satisfied for 
Hi <C i/"'' ^ 0(10 — 100 keV). This condition on Hi has non-trivial implications for the construction of explicit 
inflationary models, since values of Hi of this magnitude tend to be rather non-generic [65l| among typical classes 
of inflationary potentials. However, as discussed in Ref. [3|, such a scale for Hi is certainly not excluded (see, e.g., 
Rcfs. [la, ISJI). Moreover, a small value for Hi flts naturally within the context of the LTR cosmology. 



M. Other Astrophysical Constraints on Light Axions 



In addition to the constraints we have discussed above, there exist a number of additional astrophysical and 
cosmological bounds on theories involving light axions and axion-like particles. Indeed, particles of this sort can give 
rise to a number of potentially observable effects Q, such as a rotation of the CMB polarization, modifications of 
the matter power spectrum, and the enhanced spindown of rotating black holes. However, in order to give rise to 
observable effects of this sort, the particle in question must be exceedingly light, with a mass m < 10~^° cV. In the 
extra-dimensional scenario we are discussing here, the Newton's-law-modification constraint on the compactification 
scale Mc stated in Eq. (|4.ip implies that all a\ in the tower have masses A > 10^'^ cV in any scenario in which 
J/ < 1, i.e., in which the full tower of a\ contributes significantly to fltot- Consequently, the additional constraints on 
ultra-light axions and axion-like fields discussed in Ref. J9|] are not relevant for our bulk-axion DDM model. 
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SYNTHESIS: COMBINED PHENOMENOLOGICAL CONSTRAINTS ON AXION MODELS OF 

DYNAMICAL DARK MATTER 



In the previous section, we enumerated the individual astrophysical, phenomenological, and cosmological consider- 
ations which potentiahy constrain our bulk-axion DDM model, and we evaluated the restrictions that each placed on 
the parameter space of this model. In this section, we summarize how these individual results, taken together, serve 
to constrain that parameter space. Our particular interest concerns the preferred region of parameter space outlined 
in Ref. 0, namely fx ~ 10^^ - lO^^ Q^y, Ac ~ 10^ - 10^ GeV, and M^ chosen sufficiently small that y<l. Indeed, 
this is the region within which the full KK tower contributes non-trivially to the total dark-matter relic abundance. 

In Fig. [T21 we show the combined exclusion regions for a purely photonic axion with c^ — 1 for Aq = 1 GeV (left 
panel), Aq = 1 TeV (center panel), and Aq = 100 TeV (right panel). The shaded regions displayed in each of the plots 
are excluded by the various considerations discussed in Sect. |lVl Specifically, the exclusion regions appearing in these 
panels are those associated with helioscope limits on solar axion production (red), collider considerations (magenta). 
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FIG. 12: Exclusion regions associated with all applicable phenomenological constraints discussed in this paper for our bulk- 
axion DDM model with Aq = 1 GeV (left panel), Ag = 1 TeV (center panel), and Ag = 100 TeV (right panel). In each 
case, we have taken ^ = gc = 1, Trh = 5 MeV, and Hi = W~^ GeV, and we have assumed that the axion only couples 
to the photon field with c-, = 1. The shaded regions are respectively excluded by data from helioscope measurements (red), 
collider considerations (magenta), tests of Newton's-law modifications via Eotvos-type experiments (purple), measurements of 
the diffuse extragalactic X-ray and gamma- ray spectra (orange), observations of the lifetimes of globular-cluster stars (yellow), 
energy- loss limits from supernova SN1987A (cyan), the model-consistency requirement Aq < fx (gray), overproduction of 
thermal axions (green) , and the upper bound on the dark-matter relic abundance from WMAP (brown) . The dashed black line 
corresponds to y = tt; smaller values of y correspond to the region below and to the left of this line. 
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FIG. 13: Same as Fig. 1121 but for a "hadronic" axion — i.e., an axion which couples to both photons and gluons (and hence to 
pions, nucleons, and other hadrons), but not directly to SM quarks or leptons. For these panels, we have taken Cg — Cj = 1. 
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tests of Newton's-law modifications via Eotvos-type experiments (purple), measurements of the diffuse extragalactic 
X-ray and gamma-ray spectra (orange), observations of the hfetimes of globular-cluster stars (yellow), energy- loss 
limits from supernova SN1987A (cyan), the model-consistency requirement Aq < fx discussed in Ref. [2| (gray), and 
the 3ct upper bound on the dark-matter relic abundance from WMAP (brown). The additional requirement that the 
relic abundance be primarily determined by the misalignment mechanism (as envisioned in our DDM model) excludes 
the green-shaded region, within which a substantial population of ax is generated via interactions with SM particles 
in the thermal bath. The remaining unshaded regions of parameter space are the regions within which our DDM 
model is consistent with all of these constraints. The dashed black line indicates the contour y — n; smaller values of 
y correspond to the region below and to the left of this line. As discussed in Ref. [4I, we are particularly interested in 
the unshaded region of parameter space which falls below and to the left of this line, since this is the region within 
which not only are all of the aforementioned constraints satisfied, but also the full tower of a\ contributes non-trivially 
to r^tot- 

As we see in Fig. [121 for small Aq the most stringent constraint on the parameter space of our model is the one 
derived from energy-loss limits from SN1987A. The constraint from globular-cluster stars is also reasonably stringent, 
and the constraint derived from missing-energy processes such as pp — ?► j+I^t at the LHC is estimated to be of roughly 
the same order. However, as the y = n contour superimposed over each panel in Fig. [T^] indicates, the full tower of 
a\ contributes significantly to iltot for all Aq ^ 100 GeV. Indeed, this is precisely the Aq regime associated with the 
preferred region of parameter space for our model. We therefore conclude that within this region of parameter space, 
a photonic bulk-axion DDM ensemble constitutes a viable dark-matter candidate. 

In Fig.[Tni we consider all of the same constraints as in Fig.[T21 but for the case of a hadronic axion with Cg = c-^ = 1. 
In this case, since the a\ couple to hadrons, the constraints from SN1987A and from axion production via interactions 
among the SM particles in the radiation bath both become even more stringent. Again, as in the photonic-axion 
case, we find that the leading constraint for small A^ is that from SN1987A, and that as Aq increases, the model- 
consistency constraint becomes increasingly stringent. However, as in the photonic-axion case, we see that within the 
preferred region of parameter space for our model, a hadronic bulk axion is also consistent with experimental and 
observational limits. Thus a hadronic bulk-axion DDM ensemble is a viable dark-matter candidate as well. 

We also observe that the exclusion contours in Figs. [T^and[T51 associated with SN1987A energy-loss limits, globular- 
cluster-star evolution, collider constraints, and axion overproduction from SM particles in the radiation bath have the 
same slope. This is because all of these constraints involve the production of light axions which are never directly 
detected, and thus involve physical processes whose amplitudes include a single coupling factor between the a\ and a 
pair of SM fields. By contrast, the slopes of the constraint contours associated with other classes of physical processes 
can be quite different. The helioscope-constraint contour, for example, is related to processes in which axions are both 
produced and subsequently detected via their interactions with SM fields. Likewise, the contour associated with limits 
on features in the diffuse X-ray and gamma-ray backgrounds is due to processes involving the decays of a preexisting 
cosmological population of axions, and therefore depends not only on the couplings of the a\ to SM fields, but to their 
relative abundances as well. The slopes of these constraint contours consequently differ from those which characterize 
the contours associated with SN1987A energy-loss limits, globular-cluster-star evolution, and so forth. 

VI. DISCUSSION AND CONCLUSIONS 

In Ref. [1] , we proposed a new framework for dark-matter physics which we call "dynamical dark matter" (DDM) . 
The fundamental idea underpinning DDM is that the requirement of stability is replaced by a delicate balancing 
between lifetimes and cosmological abundances across a vast ensemble of individual dark-matter components. If 
Ref. [1[, we developed the general theoretical features of this new framework. By contrast, in Ref. [4I, we presented a 
"proof of concept," namely an explicit realization of the DDM framework in which the DDM ensemble is realized as 
the infinite tower of KK excitations of an axion-like field propagating in the bulk of large extra spacetime dimensions. 

In this paper, we have completed this study by systematically investigating all of the experimental, astrophysical, 
and cosmological constraints which apply to this DDM model. Some of these constraints pertain to theories with large 
extra dimensions in general, while others pertain specifically to our model. Among the bounds we have considered are 
constraints from limits on a\ production by astrophysical sources such as stars and supernovae; constraints related 
to the effects of late relic-axion decays on BBN, the CMB, and the diffuse X-ray and gamma-ray backgrounds; 
collider constraints on missing-energy processes such as pp — > J + I^t and pp — > 7 + I^t', constraints on isocurvature 
perturbations generated as a consequence of misalignment production; constraints on the production of relativistic 
axions due to interactions in the thermal bath after infiation; and constraints on the direct detection of dark axions 
by microwave- cavity detectors and other, similar instruments. We have verified that all of these constraints are 
satisfied within the preferred region of parameter space for our model — namely, that in which the bulk-axion 
DDM ensemble accounts for the observed dark-matter relic abundance, while at the same time the full tower of 
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axion modes contributes meaningfully to that abundance. We therefore conclude that this bulk-axion DDM model is 
indeed phenomenologically viable, and that the overall DDM framework is a self-consistent alternative to traditional 
approaches to the dark-matter problem. 

While the focus of this paper has been on the specific bulk-axion DDM model presented in Ref. 2] , we note that 
many of our results, and in many places our entire methodology, have a far wider range of applicability. For example, 
much of the formalism developed in Sect. IIVI for evaluating the cosmological constraints on decaying dark matter in 
our bulk-axion DDM model is applicable to any model in which the dark sector comprises a large number of fields. 
This is true for issues as diverse as BBN, diffuse photon backgrounds, or stellar cooling. Likewise, irrespective of 
issues pertaining to dark-matter physics, many of our results and techniques may have applicability to theories with 
large numbers of axions, such as the recently discussed "axiverse" theories [^ISal- Thus, we believe that the methods 
developed and employed in this paper can serve as a prototype for future phenomenological studies of not only the 
DDM framework, but also, more generally, any theories in which there exist large numbers of interacting and decaying 
particles. 
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